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Abstract 

We consider random Schrodinger equations on for d > 3 with a homogeneous 
Anderson-Poisson type random potential. Denote by A the coupling constant and di the 
solution with initial data do- The space and time variables scale as x 
A“^“'' with 0 < K < Ko{d). We prove that, in the limit A —> 0, the expectation of 
the Wigner distribution of dt converges weakly to the solution of a heat equation in 
the space variable x for arbitrary initial data. The proof is based on a rigorous 
analysis of Feynman diagrams. In the companion paper m the analysis of the non¬ 
repetition diagrams was presented. In this paper we complete the proof by estimating 
the recollision diagrams and showing that the main terms, i.e. the ladder diagrams with 
renormalized propagator, converge to the heat equation. 
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1 Introduction 


The Schrodinger equation is time reversible and has no dissipation. The long time 
dynamics of a quantum particle in a small random environment nevertheless exhibits 
a stochastic behavior that can be described by a diffusion equation. Our goal is to 
establish this fact rigorously in the weak coupling regime. We have announced in cni 
that under a scaling of space and time with inverse powers of the coupling constant A, 
the Wigner distribution of the solution to the random Schrodinger equation converges to 
the solution of a heat equation as A —> 0. Our approach is based on graphical expansion 
methods coupled with a certain truncation scheme. The first part of the proof was given 
in cni; the current paper contains the second and final part. To help the orientation of 
the reader, we now summarize the important notations and the main result below. 

The quantum dynamics of a single particle in a random potential is given by the 
Schrodinger equation 


= teR. (1.1) 

The Hamiltonian is a Schrodinger operator, 

H:=-^A + XV, ( 1 . 2 ) 

acting on L^(M'^) with a random potential V = Vi^{x) and a small positive coupling 

constant A. The potential is given by 

Vuix) := [ B{x - y)dfi^{y) , (1.3) 

where H is a single site potential profile and is a Poisson point process on with 
homogeneous unit density and with independent, identically distributed random coupling 
constants. More precisely, for almost all realizations co consists of a countable, locally 
finite collection of points, {y^{uj) : 7 = 1,2 ,...}, and random charges {v^{ijj) : 7 = 
1 , 2 ,...} such that the random measure is given by 

00 

= '^v^{u))5y^^^) (1.4) 

7=1 

where 5y denotes the Dirac mass at y G R'^. The Poisson process {y-^{ijj)} is independent 
of the charges {u.y('^)}- The charges are real i.i.d. random variables with distribution 
and with moments := E„ satisfying 

m2 = 1 , m2d < 00 , mi = m 3 = ms = 0 . (1-5) 
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The expectation with respect to the random process {y-y,v-y} is denoted by E. For the 
single-site potential, we assume that 5 is a spherically symmetric Schwarz function with 
0 in the support of its Fourier transform. 

A quantum wave 'ip G function can be represented on the phase space by its 

Wigner transform 

W^{x, v) := J - |)d7? , (1.6) 

with the convention that integrals without explicit domains denote integration over 
with respect to the Lebesgue measure. Define the rescaled Wigner distribution as 

W^{X,V):=e-'^W^(^^,v) . (1.7) 

The Fourier transform of the kinetic energy (dispersion relation) is e{p) := the 

velocity is given by ^Ve(p) = ^p. 

For any function /i : —> C and energy value e > 0 we introduce the notation 

[h]{e) := J h{v)6{e - e{v))dv := h{q) > (1-8) 

where du{q) is the restriction of the Lebesgue measure onto the energy surface Sg := 
{q : e{q) = e} that is the ball of radius The normalization of the measure [-Jg is 

given by 

[l](e) :=Crf_i(2e)i-i , (1.9) 

where Cd-i is the volume of the unit sphere 5"^“^. 

We consider the scaling 

t = X-^{X-'^T), x = X-^/^{X-^X) = X/e, e = (uq) 

with some k > 0. On the kinetic scale, k = 0, the limiting equation is the linear 
Boltzmann equation with a collision kernel 

a{u, v) := 27r|i3(M — v)\^5{e{u) — e(v)) . (1-11) 

The generator of the corresponding momentum jump process v(t) on the energy surface 
Sg, e > 0, is 

Lef(v} := J dua(u,v)[f(u) - f(v)], e(v) = e . (1.12) 

The choice k > 0 corresponds to the long time limit of the Boltzmann equation with dif¬ 
fusive scaling. The Markov process {'y(t)}t>o with generator Lg is exponentially mixing 
(see Lemma EH in the Appendix). Let 

^b(e) := / Tg[uW(t)i;(j)(0)]dt , v = i, j = 1,2,... d, 

Jo 
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be the velocity autocorrelation matrix, where denotes the expectation with respect 
to this Markov process in equilibrium. By the spherical symmetry of B and e{U), the 
autocorrelation matrix is constant times the identity: 

1 

Dij{e) = 5ij, := / Se[v{t) ■ v{0)]dt . (1.13) 

[ZTT) a Jq 

The main result is the following theorem. 


Theorem 1.1 Let d > 3 and i/jq G be a normalized initial wave function. Let 

V’(t) := exp{—itH)'tpQ solve the Schrodinger equation ifl.ll) . Let 0{x,v) be a Schwarz 
function on x M'^. For any e > 0, let f be the solution to the heat equation 

dTf{T,X,e) = De Axf{T,X,e) (1.14) 


with the initial condition 


fiO,X,e) :=6iX) \Mv)\^ (e) 


Then there exist 0 < Ko{d) < 2 such that for 0 < n < Ko{d) and for s and A related by 
the rescaled Wigner distribution satisfies 


lim j dX jdv OiX,v)BW^^^.,.,^-^iX,v) = ^ dX Jdv OiX,v)f{T,X,eiv)) , (1.15) 

and the limit is uniform on T £ [0, Tq] with any fixed Tq. In d = 3 one can choose 
ko{3) = 1/500. 

Our result shows that the quantum dynamics on the time scale t ~ is given 

by a heat equation and the diffusion coefficient can be computed from the long time 
behavior of the underlying Boltzmann dynamics. Heuristically, this statement can be 
understood from two facts. First, the Boltzmann equation correctly describes the limit 
of the quantum evolution on the kinetic time scale t ~ A“^, A —> 0 (see jTS], [7j, [S] 
and m)- Second, the long time limit of the linear Boltzmann evolution is a diffusion. 
This two step limiting argument is, however, misleading (e.g. in d = 2, localization is 
expected to occur at all values of A, so that no diffusion occurs). In higher dimensions, 
quantum correlations that are small on the kinetic scale and are neglected in the first 
limit may become important on the longer time scale, too. To prove Theorem II.IL we 
need to control the full quantum dynamics up to the time scale t ~ A“^“'^ and prove 
that the quantum correlations are not sufficiently strong to destroy the heuristic picture. 


The approach of this paper applies also to lattice models and yields a derivation 
of Brownian motion from the Anderson model BU- The dynamics of the Anderson 
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model was postulated by Anderson jH] to be localized for large coupling constant A and 
extended for small coupling constant (away from the band edges and in dimension d >3). 
The localization conjecture was first established rigorously by Goldsheid, Molchanov 
and Pastur m in one dimension, by Frohlich-Spencer m, and later by Aizenman- 
Molchanov pQ in several dimensions, and many other works have since contributed to 
this field. 

The progress for the extended state regime, however, has been limited. It was proved 
by Klein |14] that all eigenfuctions are extended on the Bethe lattice (see also (21 ). In 

Euclidean space, Schlag, Shubin and Wolff dZl proved that the eigenfunctions cannot be 
localized in a region smaller than for some 6 > 0 in d = 2. Chen |S], extended this 

result to all dimensions d >2 and <5 = 0 (with logarithmic corrections). Extended states 
for Schrodinger equation with a sufficiently decaying random potential were proven by 
Rodnianski and Schlag dH and Bourgain (^ (see also (HI). However, all known results 
for Anderson-type models in Euclidean space are in regions where the dynamics have 
typically finitely many effective collisions. Under the diffusive scaling of this paper, see 
(fTTOl) . the number of effective scatterings is a negative fractional power of the scaling 
parameter. In particular, it goes to infinity in the scaling limit, as it should be the case 
if we aim to obtain a Brownian motion. 

Acknowledgement. The authors are grateful for the financial support and hospital¬ 
ity of the Erwin Schrodinger Institut, Vienna, Max Planck Institut, Leipzig, Stanford 
University and Harvard University, where part of this work has been done. 

2 Summary of Part I 

2.1 Notations 

We introduce a few notations. The letters x,y,z will denote configuration space vari¬ 
ables, while p, q, r, u, v, w will be used for d-dimensional momentum variables. The norm 
II • II denotes the standard L^(M'^) norm and we set 

ll/IU,»^= E IIW’”»“/(a^)lloo 

\a\<n 

with (x) := (2 + x^Y^‘^ (here a is a multiindex). The bracket (•, •) denotes the standard 
scalar product on L^(R'^) and (•, •) will denote the pairing between the Schwarz space 
and its dual on the phase space x M'^. The Fourier transform is denoted by hat. 
For functions on the phase space, f{x,v), G M'^, Fourier transform will always be 
understood in the first variable only: 

mv) := I e-2-«-/(x,u)du. 
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By the regularization argument in Section 3.2 of m, we can cutoff the high momentum 
regime of ipQ and we can multiply B by a A-dependent cutoff function in momentum 
space at a negligible error. When dealing with estimates for any fixed A > 0, we can 
thus assume that 


suppV'o is compact, suppi3(p) C {|p| < A (2-1) 

for any fixed <5 > 0. Universal constants and constants that depend only on the dimension 
d, on the final time Tq and on V’o and B will be denoted by C. The same applies to the 
hidden constants in the O(-) and o(-) notations. 

In |inj the self-energy operator was defined as the multiplication operator in momen¬ 
tum space 


6{p) := 0(e(p)), 0(a) ;= lim 0 ^( 0 ) , Oeia) := 0£(a,r) (2.2) 

e—>0-1- 


for any r with e(r) = a, where 


0e(a,r) 


\B{q — r)pdg 
a — e{q) + ie 


(2.3) 


The function 0(a) is Holder continuous with exponent and it decays as 
(Lemma 3.1 and 3.2 in HO]). If we write 0(a) = TZ{a) — where Tl{a) and I{a) 

are real functions, and recall Im{x + i0)“^ = —tt5{x), we have 


J(a) 


—Im 0(a) = TT 


6{e{q) - a)\B{q - r)\‘^dq 


(2.4) 


for any r satisfying a = e(r). 

The dispersion relation was renormalized by adding the self-energy term: 

uj{p) := e{p) + \^0{p) , 

and the Hamiltonian was rewritten as 

H = Hq + V, Ho:=uj{p), V := XV - X^e{p) . (2.5) 

The renormalization compensates for the immediate recollisions in the Duhamel expan¬ 
sion (see Section n. The rate of the immediate recollisions is of order A^, thus their 
total effect is X^t ^ 1. The renormalization removes this instability. We note that uj{p) 
is not the self-consistently renormalized dispersion relation, but only its approximation 
up to O(A^). Since this error is negligible on our time scale due to A^t <C 1, we use uj to 
simplify the technicalities associated with the analysis of the self-consistent dispersion 
relation. 
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We define the renormalized propagator (with r/-regularization); 


Rn{a,v) 


1 

a — uj{v) + irj 


In Add endix IB ■ 1 1 we prove the following estimates on the renormalized propagator. 


Lemma 2.1 Suppose that > r] > with K < 1/12. Then we have, 


and for 0 < a < 1 


, .... ^ _ 
a — u!{p) + iri\ 

(a)V7|,|-y^) 

\h{p-q)\dp 

/ CaWhhdfl A-2(l-“) 

\a — oj{p) + fry 

-(a)“/7|,|-y2H> 

\h{p-q)\dp 

^ CaWhhdpr^^^-^^ 


|a-e(p)+fr/|2-“ (a)^/^{\q\ - ^/^\) 


( 2 . 6 ) 


(2.7) 

( 2 . 8 ) 


For a = 0 and with h := B, the following more precise estimate holds. There exists a 
constant Cq, depending only on finitely many norms, such that 


f X‘^\B{p — g)p dp 

J \a — uJ{p) — fr/P 


<l + CoX-^^^[X + \a-w{q)\^/^] . 


(2.9) 


2.2 The expansion 

We expand the unitary kernel of 77 = Hq + V (see (ESI)) by the Duhamel formula and 
after taking the expectation, we organize the expansion into sums of Feynman diagrams. 
In order to avoid the infinite summations dll in the expansion, we have reduced the 
problem to a large finite box, Ai = [—L/2, L/Tf' C with periodic boundary conditions 
(see Section 3.3 of uni). The infinite volume Poisson process was replaced with 

M 

7=1 

where M is a Poisson random variable with mean |A/^|, the points are uniformly 

distributed on A^ and the real charges u/, have distribution P^. All random variables 
are independent. Lemma 3.4 of cni guarantees that these modifications have no effect 
on the final result if L —> oo is taken before any other limit. 

After the Duhamel expansion, taking the expectation and letting L —> oo, we regain 
the infinite volume formulas for the Feynman graphs. In m we used primes to denote 
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the restricted quantities, but to avoid the heavy notation here we will omit them, except 
when stating the theorems. All quantities throughout this section are understood on 
with M random points. 

We recall the Duhamel formula from Section 4 of m- For any fixed integer A^ > 1 


N-l 


ipt ■■= e ^ V’n(i) + , 

n=0 


with 


^n(t) 

4'Ar(t) 


with the notation 


:=( 

:=( 






( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


Substituting V = —\^9{p) + with V.^{x) := VyB{x — y-^), the terms in (tmi) 

and are summations over collision history. Denote by F^, n < oo, the set of 

sequences 

7 = (7i,72,---,7n), 7j G {1,2,... ,M} U {?9} (2.13) 

and by Wly the associated potential 


:= 


XV^ 

-x^e{p) 


if 7G{1 ,...,M} 
if 7 = . 


The tilde refers to the fact that the additional I'd} symbol is also allowed. An element 
7 j of the sequence 7 G F is called potential label and j is called potential index if 7 ^ G 
{1,2,...M}, otherwise they are called id-label and d-index, respectively. A potential 
label carries a factor A, a "d-label carries }?. 

For any 7 G F„ we dehne the following fully expanded wave function with truncation 

:= (2.14) 

Jo 

and without truncation 

f [dsj]”+^ • (2.15) 

Jo 
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In the notation the star (*) will always refer to truncated functions. Each term 
along the expansion procedure is characterized by its order n and by a sequence 7 £ r„. 
The main terms are given by non-repetitive sequences that contain only potential labels, 
i.e. we define 

rr := {7 = ( 71 ,- • • ,7fc) : Ij e {!,• • ■,M}, 7 * ^ if i C . (2.16) 

The sum of the corresponding elementary wave functions is denoted by 

■■= E ^77 • (2-17) 

76rr 

The rate of collisions is O(A^), thus the total number of collisions is typically of order 
k A^t. We thus set 

K := [A“^(A2t)] (2.18) 

([ • ] denotes integer part) to be an upper threshold for the number of collisions in the 
expanded terms. Here 6 = 5{k) > 0 is a small positive number to be fixed later on. 


2.3 Structure of the proof. 

In Section 5 of m the Main Theorem was proved from three key theorems. For complete¬ 
ness, we repeat these three statements here. Recall that the prime indicates restriction 
to Al and hence dependence on L, M 


Theorem 2.2 (L^-estimate of the error terms) Let t = 0(A ^ and K given by 
(ITTKl) . If K < Ko{d) and 5 is sufficiently small (depending only on k), then 


K-l 


lim lim 
A— L^cx. 


E' 


V’t - E ^7 


k=0 


= 0 . 


In d = 3 dimensions, one can choose ko{3) = 

Theorem 2.3 (Only the ladder diagram contributes) Let n < e = 

T] = A^+'^j t = and K given by (12.181) . For a sufficiently small positive S and 

for any 1 < k < K we have 

lim = ^A(i,A:)+ (2.19) 

L —>00 ’ V / 


lim {dL,F.'Wf,.r)L = Wx{t,k,0) + 0{\^-^- 
L^oo \ 




( 2 . 20 ) 
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as A <C 1. Here 


\2kp2tr} p poo p , \ ^ 

y\{t,k) := J J dad/3 y (lld^i) \^o{pi)\‘^ 

fc+1 k 

X Yl Rvi'=^,Pj)RviPHj)'[\\B{pj+i - Pj)f (2.21) 

i=i 3=1 


\2k„2tr] j- roc j- r- K+1 ^ __— 

Wx{t,k,0) J j dad/3 y dC J (Hdt;, )0{^,Vk+i)W^^i^,vi) 

k+l - -T — ^ k 

X n R,(^a,v, + - f) n ■ (2.22) 

3=1 3=1 

The definition (IT^ does not apply literally to the free evolution term k = 0; this term 
is defined separately: 

Wx{t, k = 0,O):= J d^du d{^, v)Wo{eC, v) . (2.23) 


Theorem 2.4 (The ladder diagram converges to the heat equation) Under the 
conditions of Theorem \2JA and setting t = we have 


K-l 

lim y^\{t,k,0) 


I dX I dvO{X,v)f{T,X,e{v)) , 


(2.24) 


where f is the solution to the heat equation dm. 

The main result of [1 Oj was the proof of Theorem 12 .3L In Sections OHSl of this paper 
we prove Theorem 12.21 and in Section 1^1 we prove Theorem l2.4[ We now explain the ideas 
behind these theorems. The actual estimates are somewhat weaker than the heuristics 
predicts. 

The proof of Theorem 12.21 consists in controlling the wave functions of collision his¬ 
tories that contain ??-labels or repeated potential labels. The repeated potential labels 
correspond to recollisions with the same obstacle. Immediate recollision with the same 
obstacle occurs with an amplitude O(A^). Over the total history of the evolution, this 
would yield a large contribution of order X^t ^ 1. However, the wave functions of these 
collision histories will be resummed with those containing ?9-labels. Thanks to the choice 
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of the renormalization counterterm, X‘^6{p), the contributions of the immediate recolli¬ 
sions and the ??-labels cancel each other up to leading order. Each resummed term thus 
has an amplitude O(A^t). The full propagator in the error term, ( 12 . 121 ) . however, 

will be estimated by unitarity (|13I)- This estimate effectively loses an extra t factor. To 
compensate for it, we have to continue the expansion up to two immediate recollisions 
in the error term. 

The amplitudes of the non-immediate recollisions are much smaller and they can be 
estimated individually. Heuristically, the probability of such recollisions can be under¬ 
stood in classical mechanics. Since the mean free path is A“^, returning to an already 
visited obstacle after visiting another obstacle at distance has probability 0(A^'^). 
Another scenario is when the particle collides with obstacle 71, then it bounces back 
from a nearby obstacle 72, I71 — 72I = 0(1), and then it recollides with 71. This situ¬ 
ation is atypical and it is penalized by O(A^) because the time elapsed between these 
collisions is 0(1) while the collision cross-section is O(A^). In conclusion, the probability 
of a non-immediate recollision among the total k ~ X^t collisions is at most 0{kX‘^), thus 
the total effect of these recollisions on our time scale is negligible, even when multiplied 
with the additional factor t from the unitarity estimates for 'kAr(t). 

This outline neglects the key analytic difficulty originated from the growth of the 
combinatorics of Feynman diagrams. The amplitude of non-repetition wave functions 
can be written as a sum of fc! Feynman diagrams. Only one of them, the ladder diagram, 
contributes to the heat equation. All other diagrams can be estimated by O(A^) due to 
phase cancellations. This estimate is not sufficient to sum up all diagrams since their 
number, k\ ~ exp(A““"^*), is exponentially large. The size of a few combinatorially 
simple diagrams is indeed O(A^), but much stronger estimates were obtained in cni as 
the combinatorial complexity of the diagram increases. This improvement balances the 
increased combinatorial factor for more complicated diagrams and it allows us to control 
the expansion for non-repetition wave functions for time scale t ~ X~^~'^. 

In this paper, we extend the classification scheme to include all diagrams arising 
from the Duhamel expansion. Thanks to the stopping rules, this part involves only a 
few extra collisions. The main idea of this paper is to design a surgery of Feynman 
diagrams so that a general diagram can be decomposed into a repetition and a non¬ 
repetition part: the repetition part involves only a few variables and the integration 
can be estimated accurately; the non-repetition part is reduced to Theorem 12.111 The 
errors from the surgery are controlled by the small factors from the repetition part. This 
renders all repetition diagrams negligible. Thus we prove that among all diagrams only 
the ladder diagrams without repetition contribute to the final heat equation. In practice, 
the scheme used in this paper is much more complicated than is stated here. But this 
description gives a good first idea. 

Finally, the proof of Theorem [Tl] is a fairly explicit but delicate calculation involving 
singular integrals. The proof shows how the long time evolution of the Boltzmann 
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Figure 1: Gates and skeletons 
equation emerges from the ladder diagrams. 


3 The stopping rules 

We use the Duhamel expansion to identify the non-repetition error terms to be estimated 
in Theorem rni This method allows for the flexibility that at every new term of the 
expansion we perform the separation into elementary waves, we can decide 

whether we want to stop (keeping the full propagator as in (12.121) ') or we continue to 
expand that term further. This decision will depend on the collision history, 7, and it will 
be given by a precise algorithm, the stopping rules. The key idea is that once the collision 
history 7 is “sufficiently” atypical, i.e. it contains either atypical recollision or too many 
collisions, we stop the expansion for that elementary wave function immediately to reduce 
the number and the complexity of the expanded terms. 

Not every recollision is atypical. An immediate second collision with the same obsta¬ 
cle contributes to the main term; this is actually the reason why the dispersion relation 
had to be corrected with the self-energy X'^9{p). 

In a sequence 7 we thus identify the immediate recollisions inductively starting from 
7i (due to their graphical picture, they are also called gates). The gates must involve 
potential labels and not i ). For example, the sequence 7 = (a, id , a, b, b, c, d, d, id, d , e, e, f) 
has three gates (see Fig. P). In the sequence {a,b,b,c,c,c) there are two gates. Any 
potential label which does not belong to a gate will be called skeleton label. The index 
j of a skeleton label 7j in 7 is called skeleton index. The set of skeleton indices is S{j). 
Similar terminology is used for the gates. In the first example 1,3,6,13 are skeleton 
indices and a,a,c,f are skeleton labels, in the second example 1,6 are skeleton indices 
and a, c are skeleton labels. The d terms are never part of the skeleton. 

This definition is recursive so we can identify skeleton indices successively along the 
expansion procedure. Notice that a skeleton index may become a gate index at a later 
stage of the expansion, but never the other way around. 

The formal definition is as follows. Let In := {1, 2,... , n}. 
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Definition 3.1 (Skeleton labels and indices) Let 7 = (71,72, • • • ,7n) £ fn o,nd let 
1 * '■= (71,72, ■ ■ ■, 7n-i) its truncation. The set of skeleton indices of S { j ) C In, is 
defined inductively (on the length of as follows. If G Fi, then S{fi) := { 1 } if H 
and £'(7) ;= 0 otherwise. Furthermore, for any 7 G F^,, n> 2 , let 

( <5(7*) if % = d 

£'(7) := < S'(7*)\{n-1} if n-leS{j*) and % = fin-i 

[ S'(7*)U{n} if 7n / and [% or 

Finally, 7^ is called skeleton label ifnGS{'y). 


For any 7 G F„, let 


Hi) ■= |5'(7)I, 


be the number of skeleton indices in 7, let 


iHl) ■= {i ■ Ij = 

be the set of 0 -indices and t{fi) := |F^(7)|. Denote the total A^-power collected from 
non-skeleton indices by 

kH)] + t{fi) ( 3 . 25 ) 

Notice that r(fi) is integer. 

The exact stopping rule requires somewhat tedious definitions of different types of 
elementary wave functions. First we give these definitions intuitively, state our final 
representation formula for fit using these concepts, then we give the precise definitions 
and prove the formula. 


Sequences where the only repetitions in potential labels occur within the gates are 
called non-repetitive sequences. A special case is the set of non-repetitive sequences, 
F^^, without gates and 0 -labels. The repetitive sequences are divided into the following 
categories (Fig. | 2 I). If two non-neighboring skeleton labels coincide, then the collision 
history includes a genuine (non-immediate) recollision. If a skeleton label coincides with 
a gate label, then we have a triple collision of the same obstacle. If two neighboring 
skeleton labels coincide and are not immediate recollisions because there are gates or d’s 
in between, then we have a nest. 

We stop the expansion at an elementary truncated wave function (ItTHi characterized 
by 7, if any of the following happens (precise definitions will be given in Definition I 3 . 3 |l . 

• The number of skeleton indices in 7 reaches K (see (IrTsli l. We denote the sum of 
the truncated elementary non-repetitive wave functions up to time s with at most one 
power from the non-skeleton labels or ??’s and with K skeleton indices by 
The superscript (< 1 ) refers the number of collected A^ powers from non-skeleton labels. 
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# Skeleton index 
O Non-skeleton index 

Figure 2: Repetition patterns 

• We have collected from non-skeleton labels. We denoted by the sum 

of the truncated elementary wave functions up to time s with two power from the 
non-skeleton indices (the word last indicates that the last A power was collected at the 
last collision). 

• We observe a repeated skeleton label, i.e., a recollision or a nest. The corresponding 
wave functions are denoted by V'lf 

• We observe three identical potential labels, i.e., a triple collision. The corresponding 

tri 

wave functions are denoted by ■ 

Finally, denotes the sum of non-repetitive elementary wave functions without 

truncation (i.e. up to time t) with at most one A^ power from the non-skeleton indices or 
■d’s and with k skeleton indices. In particular, the non-repetition wave functions without 
gates and -d’s (denoted by above) contribute to this sum. For notational consistence, 

we will rename '■= explicitly indicate the number of A^-powers collected 

from gates or -d’s. 

This stopping rule gives rise to the following representation. 

Proposition 3.2 [Duhamel formula] For any K > 1 we have 

i’t = ( 3 . 26 ) 

k=0 
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K 




(2),last , (<l),rec Al),nest ,{l),tri'' 

*s,/c ' '^*s,k ' ' ^*s,k 


k=0 


Proof of Proposition l,?.^ We start with the precise dehnitions. For 7 G F^ and 
£< n we introduce the notation 77 £] := ( 71 ,... ,ji) to denote the beginning segment, 
or truncation, of length d. of the sequence 7 . 

Definition 3.3 (Sets of sequences) For 0 < r < 2, k < n we let 




n,k 


:= {7 G F„ : k{j) = k, r{j) = r} 


be the set of sequences with k skeleton indices and collected from non-skeleton indices. 
Let 

_/- ^ rh) 


n,k 


■= {7 G K,k ■ [7j = Ij' ^9]^ [\j - j'\ = 1 , j,/ 0 '5(7)]} 


be the set of non-repetitive sequences. For r = 0 we have n = k and we set F^'’ := 


p(0),nr 
k,k ■ 




be the set of non-repetitive sequences whose last element is non-skeleton. Let 

fr '= U U 7:1 

k<nr<2 

be the set of all non-repetitive sequences of length n and let 

f;:={7Gf„\fr : 7[i.n-i] G J 

be the set of sequences that are repetitive, but their proper truncations are non-repetitive. 
We let 


^(r),nr 

-,k 


:= r; n {7 G : % G S{^), 3 j <n- = jj = 7^+1} 

be the set 0 /triple-collision sequences, i.e. sequences whose last entry 7 ^ is a part 
of a triple collision with a gate. Let 

■■= nn{7 G fW : 3j G 5 ( 7 ),i < n- 2 , 7 ,- = 7 ,, 5(7)n{i+l,... ,n-l} / 0} 

be the set 0 /recollision sequences. Finally, let 


W) 


\{r)^nest _ 


:= nnUeflri \(r 


S’'“u77") : 3i€S(7)J<«-2.7j=7„} 


n,k ■ n 

be the set of nested sequences. In all cases we introduce the notation 

R 

lift 


f (<R),# ._ 

^ n,k 


,,k 


r=0 


where ff = tri,rec, nest, nr, last refers to the structure of the wave function. 
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Notice that the last index n of any 7 E f * is a skeleton index, in particular H 

r* = 0. This index can create a repetition in three different ways: triple collision, 
recollision or nest. It is therefore clear from the definition, that the sets 

and for 0 < r < 2 are disjoint. Moreover, for triple collision and nested 

sequences we have r > 1. The next lemma shows that these sets include the appropriate 
beginning segment of any infinite sequence. 


Lemma 3.4 Let a positive integer K be given. Let 7 = ( 71 , 72 ) • • •) £ Loo be an infinite 
sequence. Then there exist a unique k < K and n E [/c, /c + 4] such that the truncation 
of length n of^, 77 ,n] belongs to the (disjoint) union 


f{n,K) :=f 


(<l),nr 

n,K 


u 


K 

u 

k=0 


p(<l),rec , , -p,(l),riesi 


n,k 


n,k 


ur, 


n,k 


ur 


(2),last 
n,k 


(3.27) 


Proof. We look at the increasing family of truncated sequences 7 [i, 2 ], 7 [i, 3 ])-- - in¬ 
ductively. If 7[i,n] £ some n and r < 1, then it falls into the first set of 

(ITTTl) . 

Otherwise there is an n < -t- 4 such that 7 [i^n-i] is non-repetitive with r < 1 , but 

7 [i^„] is either repetitive or ^(yji „]) = 2. If it is repetitive, then 7 ^ is a skeleton index, so 
^(7[i,n]) = ^(7[i,n-i]) r 1 , and the repetition can be a triple collision, recollision or nest 
with k < K. If 7[i,ri] is non-repetitive, then the r has increased from ^(77 = 1 to 

r( 7 [i „]) = 2 and 7 ^ is non-skeleton. These four possibilities correspond to the remaining 
sets in The disjointness of these sets follows from their definition. □ 


For 0 < r < 2 and ff = rec, nest, tri, last, nr, let 




M,# 


k+2r 

n=fc+r-gpCr),# 


be the wave function with k skeleton labels, with total power collected from non¬ 
skeleton terms and with recollision, nest, triple collision or no repetition (with the last 
collision being skeleton or not) specified by ff. The notation (*) indicates that the same 
definition is used for the wave functions with or without truncation. Finally we set 








+ F 


( 1 ),# 

{*)t,k 


We stop the expansion at the elementary truncated wave function (imi) characterized 
by 7 £ r„, if 7 falls into one of the sets in (THTl) . but none of its proper truncations 
7 [i ,i/] fell into the appropriate sets (|3.27|) with n replaced with n'. Lemma f3.4l shows 
that the expansion is stopped for every term for a unique reason. This procedure proves 
Proposition 13.21 □ 
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4 Error terms 


The content of Theorem o is that the main contribution to the wave function in 
comes from the fully expanded non-recollision terms with r = 0, i.e V’l Here 
we show that the contribution of all other terms are negligible. Each error term in 
(HT:^ has a specific reason to be small. The result can be summarized in the following 
Theorem which is proven in Sections 0] and [21 We recall that prime indicates restriction 
to Al. 


Theorem 4.1 We assume t = X ^ ^T, T e [0, Tq], and 1 < k < K. If k < 34 ^^ 39 , 
thcTi 

lim = o(A^+2''+2'5), a ^ 0, (4.1) 


L—>co 


for the following choices of parameters: = rec, r = 0 , 1 }, = nest,tri, r = 1} or 

{# = last, r = 2}. Furthermore, for k = K and r = 0,1, 


hm = o(A"+2.+25^ 


L —>00 


and for k < K, 


lim = o(A 2 ^+ 2 ^) 


r M yt^k 

L—^oc ’ 


(4.2) 


(4.3) 


Proof of Theorem, \2. ‘A using Theorem \4-l\ By (13.261) and the unitarity of the operator 
^ we have 


E' 


f- 


k<K 


K 




k=0 


^ *s,k 


(4.4) 


where the value of r is determined by ^ according to the terms on the right hand side of 
(IT^ . The non-repetition terms with r = 1 (first term on the right hand side of (I 06 l) l 
are fully expanded and there is no need for unitarity. After a Schwarz inequality, 


e 1 


K-l 

Ed 

k=0 


/ (l),nr 
k 


K-l 


< 




A"’”'I 


k=0 


Given Theorem sm all these error terms are negligible if we first take L ^ 00 and then 
A ^ 0. □ 

The natural way to prove Theorem 14.11 would be to rewrite E||^|p into a sum of 
Feynman graphs, similarly to Proposition 7.2 in m, then to identify a repetition sub¬ 
graph of a few vertices (recollisions, nests etc.) that renders the graphs small and remove 
them by graph surgeries after extracting a small factor. We then should sum up the re¬ 
maining graphs on the core indices for all possible permutations and lumps as in Section 
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9 of [ini- For the sake of technical simplifications, however, at the price of a smaller k 
we follow a somewhat different path. We first identify the vertices in the graph (called 
core indices) that carry no complication whatsoever (no repetition, no gate, no 9). We 
then symmetrize the non-core indices in and ^p,hy a, Schwarz inequality to reduce the 
number of repetition patterns. For graphs with sufficient high combinatorial complexity 
(with large joint degree, see Definition 14.41 below) we simply neglect the possible gain 
from the repetition patterns by removing them from the graph with a crude estimate. 
The necessary small factor will come from Proposition 9.2 of m with q being large. 
For graphs with low complexity, the gain comes from analyzing the repetition patterns 
case by case. Since there are not too many low-complexity graphs, we can neglect the 
possible gain from the complexity and still sum up for all combinatorial patterns of the 
core indices. 

4.1 Feynman graphs and their values 

In this section we collect the necessary definitions from [1 Dj to estimate the values of 
Feynman graphs. More details can be found in Sections 7.1 and 7.2 of m- 

The Feynman graph is an oriented circle graph on V > 2 vertices and with two 
distinguished vertices, denoted by 0, 0*. The number of vertices between 0 and 0* are 
n and n', in particular N = n + n' + 2 . The vertex set can thus be identified with the 
set V = Vn,n' ’■= {0, 1 , 2 ,... , re, 0*, re', n' — 1,..., 1 } equipped with the circular ordering. 
We set In := {1, 2,... re} and In' '■= {1, 2,... ,re'}. The set of oriented edges, T(V), can 
be partitioned into >C(V) = £ U £ so that £ contains the edges between U {0,0*} and 
£ contains the edges between /„/ U {0,0*}. 

For u G V we use the notation v — 1 and v + 1 for the vertex right before and after 

V in the circular ordering. We also denote e^_ = (u — l,u) and e„+ = (u,u -|- 1) the 
edge right before and after the vertex v, respectively. For each e G £(V), we introduce 
a momentum We and a real number Oe associated to this edge. The collection of all 
momenta is denoted by w = {we '■ e G £(V)} and dw = ^gdwe is the Lebesgue 
measure. The notation u ~ e will indicate that an edge e is adjacent to a vertex v. 

Let P = {P^ : /U G /} be a partition of the set V \ ( 0 , 0 *} = In U In' into nonempty, 
pairwise disjoint sets, where I = I(P) is the index set to label the sets in the partition. 
Let m(P) := |L(P)|. The sets are called P-lumps or just lumps. We assign an 
auxiliary variable, re^ G p, G 7(P), to each lump. The vector of auxiliary momenta 
is denoted by u := {re^ : p G 7(P)}. We will always assume that they add up to zero 

^ re^ = 0 ( 4 . 5 ) 

and that they satisfy |re^| < 0(A“^''“'^^). The set of all partitions of the vertex set 

V \ { 0 , 0 *} is denoted by Pv When we wish to indicate the re, re' dependence and 
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identify V \ {0, 0* } with U then the set of all partitions on U /„/ will be denoted 
by 'Pn,n’ instead of Vv- 
For any P C V, we let 

L+{P) := {(z;,7; + l) e £(V) : v + l^P,v€P} 

denote the set of edges that go out of P, with respect to the orientation of the circle 
graph, and similarly L^{P) denote the set of edges that go into P. We set L{P) := 
L+(P)UL_(P). 

For any ^ G we define the following product of delta functions 

A(P,w,u) := =tu;e) JI X] ±We - Uf.'j . (4.6) 

eeL±({0*}) M6/(P) e&L±(P^) 


The sign ± indicates that momenta We is added or subtracted depending whether the 
edge e is outgoing or incoming, respectively. 

For each subset ^ C V \ {0,0*}, we dehne 


Mg{-w) := \^o{We)\ \B{We,_ - We,+ )\ JJ “ ^e„+) ■ (4.7) 

e~o i;ev\{o,o*}\e v&g 

We also define the restricted Lebesgue measure 

:= l(|tc| < C)dtc , C := d/r(w) := (8)ed/x(rce) . (4.8) 

On the support of A this restriction will not substantially influence our integrals (see 
(7.9)-(7.10) of uni). 

With these notations, we define, for any P G Py and s' = 0,1,2,..., the P-value of 
the partition 


P„(P,u,q;) := A 


N-2 


sup 

g : \g\<g. 


d/r(w 


n 

ee£(V) 


1 


Oe - w(We) + ir]\ 


A(P, w,u)A/e(w) 


(4.9) 

The P-value depends also on the parameters A, r/, but we will not specify them in the 
notation. We will also need a truncated version of this definition; 


P*g P,u,q: := A^ sup /d/r w TT |- - ———— A P, w, u)A7g w . 

Q-.\Q\<9J ,tc\v) \ae-Uj{We)+ir]\ 

egr({0*}) 

(4.10) 

Let P G Vn,n' be a partition on the set In U In' ■ The lumps of a partition containing 
only one vertex will be called single lumps. The vertices 0 and 0* will not be considered 
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single lumps. Let G = G(P) be the set of edges that go into a single lump and let 
^(P) := |G(P)| be its cardinality. In case of n = n', we will use the shorter notation 
Vn = Vn,ni 'Pn = Pn,n etc. We will always have 

|n — n\ < g(P) < 4, n, n < K. (4-11) 


We also introduce a function Q that will represent the momentum dependence of the 
observable. We can assume, for convenience, that ||Q||oo < 1- We define 


M(w) := [-e{We)] [-6{We)] B{We-We+l) B{We - We+l) 

e^CDG cGCDG 


eGC\G 

e/0* 


e€C\G 

e'/'O 


X ^j;{Weo+)'^|^iWeo_)Q -(W^eo*- +^^0* + ) 


where e + 1 denotes the edge succeeding e G >C(V) in the circular ordering. 
Let a, (3 £ E., P £ Vn n', and define 


V(P,a,(3) J d/i(w) J] 


1 


tM 


1 


a - U}{We) - if] - Uj{We) + IT] 

een ee£ 


(4.12) 


(4.13) 


xA(P,w,u = 0)A4(w) . 

The truncated version, 14(P, a, /3), is defined analogously but the a and /3 denominators 
that correspond to e G L({0*}) are removed. 

We set y ;= and define 

^m(P) := ^ J J ^dadp a,/3) (4.14) 

and 

2tri r r-Y 

-E^(*)g(P,u) := J J ^dadp E^^)g{P,u,a) , (4.15) 

where a in Li(^)g(P, u, q) is dehned as ag = a for e G £ and a^. := /3 for e G C. The 
notation (*) indicates the same formulas with and without truncation. We will call these 
numbers the V-value and E-value of the partition P, or sometimes, of the corresponding 
Feynman graph. Strictly speaking, the V- and E-values depend on ^ through A and 
the P-value depends on the choice of Q as well. When necessary, we will make these 
dependencies explicit in the notation, e.g. E^ or V^(P;Q). The £-value is a convenient 
estimate for the P-value of the graph (see (7.14) of [TH] ! 

|P(,)(P)| < (CA)^L;(,),(P,u = 0) (4.16) 
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with g = (7(P). We will use the notation := u = 0). 

For the graphical representation of the Duhamel expansion we will really need 

^W(P) ^ |j^dad/3P„(P,a,/3) , (4.17) 

i.e. a version of V(*)(P) with unrestricted dad/3 integrations. (The circle superscript 
in V° will refer to the unrestricted version of P). The difference between the restricted 
and unrestricted P-values are negligible even when we sum them up for all partitions 
(Lemma 7.1 of pfU]!. 

Sometimes we will use the numerical labellingjrf the edges. We will label the edge 
between (j — l,j) by ej, the edge between {j,j — 1 ) by ej and we set Cn+i ■= (n, 0 *), 
= (0)1) ■= (1)0)- Therefore the edge set C = >C(Vn,n') is 

identified with the index set In+i U In'+i and we set 

Pj-.= We^, Pj:=We~.. (4.18) 

4.2 Resummation for core indices 

We need to identify the non-repetitive potential labels in a sequence. 

Definition 4.2 (Core of a sequence) Let 7 £ then the set of core indices of 7 

is defined as 

;= |j e S'(7) ; 7^. ^ 7., Vi 7^ j| 

and we set 0 ( 7 ) = |/“'’®( 7 )|. The corresponding 7 j labels are called core labels. The 
subsequence of core labels form an element in T”'', i.e. a sequence of different potential 
labels. The elements of 

C(7) :=/„\|/r'(7)U7"(7)l 

are called non-core potential indices. 

In other words, the core indices are those skeleton indices fPehnition l,V. /I) that do 
not participate in any recollision, gate, triple collision or nest. Given the stopping rules 
(Section|3I), the number of non-core potential indices and 0 -indices together is at most 4. 
The number of core indices c = 0 ( 7 ) is related to the number of skeleton indices k = k{f) 
as follows 

{ k if # = nr, last 

k — 1 if # = triple (4.19) 

k — 2 if ff = nest, rec. 

For any 7 £ F^, the index set In is partitioned as into core indices, 

non-core potential indices and 0-indices. Let r = r( 7 ) := (ri, r 2 ,..., Tc) £ F”^ denote the 
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core labels of the sequence 7 . We also introduce the notation T[afi] •= {'Ta,Ta+i,... ,Tb) 
if a < b, and T^a,b] = 0 otherwise. 

Now let 7 e r{n,K) U (see (Id .2 71) and Definition ESI for the notation). We 

recall that the total number of gates and 0’s is given by r. The possible values of r are 
determined by # = rec, nest, triple, last according to (lT?7l) or r = 1 if = nr and 
k < K. We will refer to a gate or 6 index as a gate/0-index in short. 

We rewrite each error term in by first summing over core labels r. For fixed 

number of core indices c we sum over all possible locations of non-core indices. If a 
non-core index is inserted between the {w — l)-th and w-ih. core indices, we characterize 
its location by w. 

The locations of non-core indices within the sequence are given by a location code 
w. For example, if # = last, then w G Ic+i encodes that the first gate/0-index is located 
between the {w — l)-th and tc-th core indices. The location of the second gate/0-index 
need not be encoded because it is fixed to be after the last core index. If ^ = rec and 
r = 0, then w € Ic encodes that the first recollision label is between the {w — l)-th 
and rc-th core indices. The most complicated case is = rec, r = 1, when the code w 
consists of two numbers, w = {wi,W2) G IcX Ic+i, where wi and W2 encode the location 
of the recollision and gate/0, respectively. If wi = W 2 , an extra binary code determines 
whether the gate/0 is immediately before or after the recollision. The set of possible 
location codes therefore depends on c and r, and it will be denoted by VF = . 

The detailed description of the set W in the other cases is obvious but lengthy and 
we omit the formal details. The precise structure of W is not important, but we remark 
that its cardinality satisfies |IF| < (c - 1 - 1)^ in all cases. We thus have the following 
resummation formula: 



Tgrj?^ waw 


where is the wave function with core labels r, with structure with r gates/ 0 - 

indices and with location of non-core indices given by w. We note that the wave function 
'll:', includes a summation over the non-core labels with the restriction that they are 
distinct from the core labels r. 


Having specified the locations of the r gates/0-indices within the sequence of core 
indices, we introduce another code h G {g,9Y) called gate-code, to specify whether 
there is a gate or a 0 at the given location. This gives the decomposition 



(4.21) 


h&{e,gY 


I h 4^ 

with the obvious definition of 'll), 
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4.3 Symmetrization of the non-core indices 

Starting from (lOl . we can use the Schwarz inequality 


It = 


E E'( E '‘'(St. E 


! h),# 


w^w'^W t'eF^^ 

/ (r),# ,/ (r),# 


loeVK T,r'er!?’' 


where c is given by k and ^ according to (ICTHl . and recall that W depends on 7 ^, c, r. 

Notice that any non-core potential label appears in pair (in a gate, nest or recollision) 
and none of them appear more than six times by the stopping rules. Since the first, 
third and fifth moments of the random variables are zero, the expectation in (ii:^ 
is nonzero only if t' and r are paired, i.e. if there is a permutation ct G ©c such that 
t' = meaning r' = Therefore 


E'liv'iSf II" < iM'i E E E E'(</>;S;t 


w^W c7E©c tEF!?'’^ 


siiriEE E E E'(<if„,-/-;SJH,„ 

w&W (Tg6c h,h'£{g,d}^ rSEJ?’’ 


(4.22) 


Note that each wave function been further decomposed into a sum over 

h-codes according to (lOTl) . However, we did not estimate the h X h' cross terms by 
an additional Schwarz inequality. The term with a gate must cancel the term with a 9 
exactly at the same location, i.e. and would not individually be negligible, but 
their sum is of smaller order. 

Notice also that each wave function ibi may involve summations over one or two 
further non-core potential labels. We use the convention that the recollision or nest label 
is denoted by v, the label of the gate or triple is denoted by /r. In case of a second gate, 
^ = last, its label will be fi. According to the non-repetition rules, {v, /r, /i) may not 
coincide with each other or with any element of r. In the products '‘p' w 
there is no repetition among r and cr(r) indices other than the ones prescribed by a. 
However, if the additional non-core labels within w denoted by ,yi' or pi', 

then there may be a few coincidences between primed and non-primed non-core labels. 
Those coincidences are allowed that do not violate the non-repetition rules requiring that 
ly, fi, jl are distinct and their primed counterparts are also distinct among themselves. 

Once the number of core indices c, a location-code w and a gate-code h are fixed, this 
defines a unique insertion of the non-core indices into the sequence of core indices Ic- The 
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core and non-core indices in the given order can be identified with where n, the total 
number of collisions, is given hy n = k + r+\{j : hj = g}\ and k is given by (14.1911 . This 
naturally defines an embedding map s = : Ic ^ In- Similarly, 's^ : 1^ In can be 

defined. The precise definition depends on 7 ^, r and h in a natural way. For illustration, 
we describe the most complicated case; # = rec, r = 1 and h = g. In this case n = c + 4 . 
For definitiveness, we consider the case when the recollision precedes the gate wi < W 2 - 
In this case the complete collision sequence is A*)^)- Let 

{ j if j < wi 

j + 1 if Wi <j < W2 
j + a if W 2 < j < c . 

With this notation, the pairing of core indices, originally determined by cr G ©c, is given 
by the pairs {s^(j), (<^0))} as subsets of the full index set /„ U In'- 

Given (^, c, a, w, h, h'), we define the partition Dq = Do(#, c, a, w, h, h') of In U In' 
by lumping only those indices that are required to carry the same potential label by the 
prescribed structure ^ and the gates. The vertices with 9 always remain a single lump, 
the remaining vertices are paired. 

For example, if # = rec, r = 1, w = {wi,W 2 ), h = h' = g, we obtain 

Do := ie/c) {^ 1 ) {^2 + 1, u >2 + 2}, {u>i, n}, {w 2 + 1 , u )2 + 2}| 

(see Fig. El and all other cases are similar. The elements of Dq consisting of pairs of 
core indices, {s{j),'s{a{j))}j^i^, are called core elements of Dq and those elements of Dq 
that contain non-core indices will be called non-core elements. In the example above, 
the last four elements are the non-core elements of Dq describing the two gates and two 
recollisions. The non-core potential labels u, g, v', g' will correspond to these non-core 
elements, respectively. 

Some of the non-core elements may be lumped together according to the possible 
coincidence between the {i',g} and {v',g''\ since the non-repetition rule do not prevent 
it. This procedure defines new partitions D that we will call derived partitions, denoted 
byD ^ Dq. In this particular case, there are seven possible lumpings of the four non-core 
elements without violating the non-repetition rule within the sets (r, u, g) and (r, v',g'). 
Figure 10 shows the partition D derived from the above partition Dq when g = g' but 
1 / / v'. 

In general D is defined by lumping together a few non-core elements of Dq under the 
constraint of the non-repetition rule (Fig. |1|. The single elements of Dq that correspond 
to 9 are never lumped. Note that each of these pairings gives rise to the appearance 
of exactly four or six identical potential labels; higher moments do not appear. The 
number of such quartets and sextets is denoted by ^ 4 (D) and ^ 6 (D). Clearly £» 4 (D) < 2 
and ^> 6 (D) < 1 . 
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Figure 3: Symmetrized recollision with a gate. Dq consists of the paired vertices 



Figure 4: Partition D lumps some non-core elements of Dq 
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Let D* C D denote the collection of non-single elements of D. Note that for each 
element of D* one selects a distinct potential label. The quantity (14.221) contains a 
summation over all such potential labels. We will use the connected graph formula 
(Lemma 6.1 from |1U| 1 for the index set D*. 

Let A G A(D*) be a partition of the set D*. We define P(A,D) G 'Pn,n' to be the 
partition of U In' whose lumps are given by the equivalence relation that two elements 
of Inkin' are P(A, D)-equivalent if their D-lump(s) are A-equivalent. The single lumps 
of D remain single in P (these are the 0 indices). 

We recall the definition of V(*)(P) and V’^°^(P) from (|4.14j) and (14.1711 . Since we will 
compute the L^-norm, the momentum shift at 0* is chosen to be ^ = 0 (see JSI)) and 
the Q function in inni), representing the observable, will he Q = 1. Furthermore, when 
defining Feynman graphs, we will always assume the following range of parameters (see 
(7.25) of ^01) unless stated otherwise: 

Tj = a2+^, t = T G [0, To], K = [A-'^(A2t)], k < K, C = g<8 

(4.23) 

with a sufficiently small <5 > 0 that is independent of A but depends on k. We recall 
that rj is the regularization of the propagator, K is the upper treshold for the number 
of skeleton indices, k, in the expansion, ( is the momentum cutoff (see (gSI) and g is 
the number of exceptional vertices where the standard \B{win — Wout)\ potential decay is 
not present (this happens for the single lumps). All estimates will be uniform in ^ and 
in T G [0,ro]. 

For each fixed {#,c,a,w,h,h'), by using (11.51) and the connected graph formula, 
similarly to Proposition 7.2 of m we obtain 


Imi ^ E'(^ 


rer" 


ih,# 




D:^Do 


^e(D) 


c(A)P(:)(P(A,D)) (4.24) 

Ae.A(D*) 


with := (we recall and (11.51) 1. The summary of the 

results in this section is 


Proposition 4.3 Let k < K, let = rec, nest, triple, last and r he one of the possible 
values allowed by or r = 1 if = nr, k < K. Let c he given by and let 

W = be the set of location codes. Then 

E'llAA’f 11^ < |IE| ^ ^ c{A)V',:,(P(A,D)) , 

wGW (tG&c h,h'G{g,0}'' D^Do AsA(D*) 

(4.25) 

□ 
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4.4 Splitting into high and low complexity regimes 

Given {^,c,a,w,h,h'), we consider the partition Dq of In U In' as defined above and 
let D Dq. Note that the collection D* contains all core elements of Dq, i.e. all pairs 
of core indices {s(j), The restriction of a partition A G A(D*) onto these 

core elements can therefore be naturally identified with a partition of Ic using the map 
{s(j), s((T(j))}jg/^ j G Ic- We denote this restricted partition by A. In the sequel we 
shall therefore view A G Ac, i.e. as a partition on Ic- 

The restricted partition A together with a also generates a partition P(A,cj) on the 
set Ic U Ic- The lumps of P(A, a) are given by U ct(A^), where are the lumps 
of A. Notice that the (s,S)-image of the restriction of P(A,D) G Vn,n' onto the set of 
core indices U is exactly P(A,(t). Since the cardinality of non-core elements 
of D* is at most 4, for any given D and A there can exist at most (c -|- 4)^ partitions, 
A G A(D*), whose restriction onto the core elements is A. 

We recall the definition of joint degree from m-- 

Definition 4.4 (i) Let A G Ak he a partition of Ik = {1,2..., A:}. Set := \Ay\, 
V G 4(A), to be the size of the v-th lump. Let 

S{A):= U A, 

v£I(A.) 

0,1/ A 2 

he the union of nontrivial lumps- The cardinality of this set, s(A) := |S'(A)|, is called 

the degree of the partition A. 

(ii) Let A G Ak and a G &k- The number 

q{A, a) ■-= max {deg(a), ^■s(A)| (4.26) 

is called the joint degree of the pair (u. A) of the permutation a and partition A. 

The sum will be split into two parts and estimated differently. In the regime 

of high combinatorial complexity, i.e., when the joint degree q{A,a) of a and A is bigger 
than a threshold q > 1 (to be determined later), then we can use the method of Section 
9 (especially Proposition 9.2) from jlOj robustly. This will be explained in Section [4.51 
For low combinatorial complexity we use the special structure given by the recollisions, 
nests, triple collisions or gates ISection 14.6B . The threshold q will be chosen differently 
for the estimates and for (ESI). 

The precise estimate is the following 

hm E'll^' yf IP < (!) + {II) + 0(A5) (4.27) 

L—>oo ’ 
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with 


{I):=\W\ic + 4fY.T.Yl sup{|y(,)(P(A,D))c(A)| : A = A'} 

w£W h,h' crG&c DJ^Dq A'eAc ^ 

q(A',a)>q 

(4.28) 

where the supremum is over all possible A £ A(D*) whose restriction A is the given 
partition A'; and 


(V) l»r| E E 

w£W cre©c 


E E E 1^(*)(P(A,D))c(A) 

h^h'^{g,6}'^ D>-Do Ag^(d*) 

q(A,(T)<q 


(4.29) 


We recall that Dq is determined by (^, c, a, w, h, h'). The error term O(A^) comes from 
replacing V(°)(- • •) with V(*)(- • •); see Lemma 7.1 of pITIj . 


4.5 Case of high combinatorial complexity 

Here we estimate the term (I) in (I4.28|) . Clearly < (m 4 )^(m 6 ) < C. We estimate 

V(*)(P(A,D)) by using (14.161) . Then, by applying Operation I from Appendix o we 
break up all the lumps in the partition P(A, D) that involve elements from non-core 
indices, U , in such a way that all non-core indices must form single lumps. Let 
P*(A,D) denote this new partition. Note that the projection of P*(A, D) onto the core 
indices is unchanged. The number of application of Operation I is at most 6. Using 
Lemma mu we can estimate i7(^,)g(P(A, D)) in terms of supuill(^,)g (P*(A,D),u) with 
an additional factor of at most A®, where 

A:= [CK(f = . 

Then we apply Operation II (Appendix 0 to remove all single lumps with non¬ 
core indices and use (mu) from Lemma mil After removing the non-core indices, the 
remaining vertex set can naturally be identified with IcU/c by using the (s, s) maps, and 
the partition P*(A, D) restricted to core indices 7“'’® U is identified with P(A, cr). 
Lemma mu is applied at most 8 times, therefore we obtain that for any a, D and A' in 
the sum (jUSHI): 


|U(,)(P(A,D))| < CK\\r^-^f sup E(^,)g{A\a,vi) . (4.30) 

u,5<8 

The application of Operation II is schematically shown on Fig. [3 

The summations over h, h' G {g, 0}^ and D Dq in H4.28I) contribute with at most a 
constant factor since r < 2 and the number of different D’s is at most 7. The cardinality 
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Figure 5: Estimate after removing all non-core indices 


of W can be bounded by (c -|- 1)^ and c< K < CX 


— K—S 


Therefore we obtain 


(I) < cA-8-d6+2d)-o{5) ^ ^ i^(^^^(A',a,u)|c(A' 


u,5<8 


(tG(Sc 

q(A',a)>q 

Using Proposition 9.2 from m, we have 

(I) < (7A'^[5“(¥'^+¥)'‘“0((5)]-8-(16+2d)K ^ 

We immediately see, that the contribution of (I) to the error term in (I4.dl) in Theorem 
EH satisfies the announced bound with a sufficiently small 6 if 


K < 


2q-48 


(34d -h 39)q + 112 -h 12d ’ 

The bounds EU-dilSl) are satisfied if 

2q-72 

^ (34d -h 39)q + 112 -h I2d ’ 

4.6 Case of small combinatorial complexity 

Here we control the term (II) in (14.291) . First we estimate the combinatorics. 
Lemma 4.5 For any q £ N, c < K and structure type we have 
sup V sup V |c(A)| < (CgA)39+3 , 

q{A,a)<q 

where we recall that Dq depends on {#,c,a,w,h,h'). 

Proof. The bound 

#{a G : ^(a) = i} < 


(4.31) 


(4.32) 


(4.33) 
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(see (8.14) from [12] ) shows that the number of permutations cr G ©c with deg(cj) < q is 
bounded by {CKY using c < K. The number of A’s whose restriction yields the same 
A is at most (c + 3)^ < (CK)'^. The number of A G Ac with s(A) < 2q is bounded by 
^ 2 q-i < Finally, |c(A)| < < i2q)^'^-^ □. 

The individual terms in are estimated in the following Proposition whose proof 

will be given in Section (Hj 

Proposition 4.6 We assume (j4.23l) and k < 34^^39 • Let a G ©c, w G h,h' G 

{g,0Y, where ^ and r vary in the different eases and let D Dq(#, c, cr, tc, h, h'). 

1) Let A G A(D*) such that q{A, a) < q, where q is a fixed number. Then the following 
individual estimates hold. 

(la) [Many collisions] Let [[ = nr, r = 0,1 and c = K, then 


|K(P(A,D))| < C-^As^ . 

(4.34) 

(lb) [Recollision]. Let # = rec, r = 0,1, then 


|K(P(A,D))| < _ 

(4.35) 

(Ic) ]Triple eollision] Let # = triple, r = 1, then 


|K(P(A,D))| < _ 

(4.36) 


2) Now let A' G Ac be given. Then the following estimates hold: 
(2a) [Non-repetition with a gate] Let [[ = nr, r = 1, then 


sup 


^ ^ ^ P(P(A,D))c(A) < . (4.37) 


D)-Do A6 ^(d*) 
A=A' 


(2b) [Last] Let [[ = last, r = 2, then 


sup 

(T.W 


^ ^ ^ 14(P(A,D))c(A) < 


h,h'&{g,eY'D)~T)o AeA(D*) 
A=A' 


(4.38) 


(2c) [Nest] Let [[ = nest, r = 1, then 


sup 

(7,W 


^ ^ ^ K(P(A,D))c(A) < . 


h,/i'e{g,0} D^Do AeA(D*) 
A=A' 


(4.39) 
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Combining Lemma 14.51 with these estimates, and using \ W\ < K^, we see that 


{II) < 

for the case # = nr, r = 1 (case (2a) above), and 

{II) < {CqK)^‘^~^'^ 

for all other cases (with q > 2). So the contributions of the error terms from (II) to 
(see (11123)) satisfy the bound (jl3I) if 


and they satisfy (|1I3-(I1I3 if 


K < 


I 

9q + I7d + 51 ’ 


(4.40) 


2 

{4d + 3)q + {12d + 9) 


(4.41) 


and 5 is sufficiently small. Combining this with (IOT]) - (l021l and optimizing, we obtain 
that there exists KQ{d) > 0 such that for any k < kq, the systems of inequalities (I4.i-ill) - 
(lOni) and have solutions for q. For d = 3, the optimal Ko{d) is a bit above 

This finishes the proof of Theorem 14. IL □. 


5 Proof of Proposition 14.6 

In each case except (la), the corresponding Feynman graph has a specific subgraph of 
a few vertices (recollision, nest etc.) that renders the value small. We shall prove that 
this subgraph gains at least a factor ^4+2re+2(5 Theorem 14.11 Then we remove 

all repetition patterns from the graph, we use the robust bounds 

sup sup£'((T, u) < CIlogAI^ 

o-eSfc u 

sup sup {a, u) < CA^IlogAp 

o-eSfe u 

from Lemma 10.2. of m to conclude the estimate. 

5.1 Many collisions 

The estimate in case (la) will come from the fact that any graph can be robustly esti¬ 
mated by the ladder graph and the value of the ladder of length L always carries a factor 
1/L!. This effect is the best seen in the time integral form. We hrst change 14 (P) back 


(5.1) 

(5.2) 
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to V^°(P) with an error smaller than by Lemma 7.1 from |1()| . We then 

apply the iL-identity (formula (6.2) in jlOj i to the definition of V)°(P) given in (14.171) to 
obtain 


F;(P) = A"+”'+^ 


dpdp K{t, p, n)K{t, p, n^)A(P, w, u = O)Al(w) 


with P = P(A,D) and 


K{t,p,n) 



We recall that g = ( 7 (P(A,D)) denotes the number of single lumps, or, equivalently, the 
number of 6 labels in h and h'. Note that we use the labelling w and p, p in parallel, 
keeping in mind the relabelling convention from (14.181) (see also Section 7.2 in |in| for 
more details). 

We use a Schwarz inequality: 


K°(P)I / / dpdp [|A(t,p,n)|2 + |A(t,p,n')PlA(P,w,u = 0)|A1(w)| 


By using Operation I, we can break up A into single lumps. Since s(A) < 2q, we have 
s(A) < 2q + 4, thus Operation I will be applied at most 2q + 3 times. 

If r = 0, i.e. the original graph was a non-repetition graph, and thus n = n' = k, 
then the trivial partition Aq corresponds to a partition P with a complete pairing. Thus 
both p and p momenta can be used as independent variables and 

k 

The estimate (UMI is then completed by the bound from the following Lemma 
with any 1/2 < a < 1. The proof will be given below. 

Lemma 5.1 For any 0 < a < 1, t = there exists a eonstant Ca sueh that 

^ (F 7^\~2—Ka\fc—1 

dp\K{t,p,k)\‘^\ipo{pi)\^Yl\B{pj -Pj+i)\‘^ < -llogAp . 

In partieular, 

I{k) < (CaA-2+'5“)^-^ (5.3) 

if k > and A <C 1. 


|K°(P)| < A2''+3(CA)2^+3 J 
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Figure 6 : Removal of gates from a recollision 

Proof. This lemma is essentially the same as Lemma 3.1 in jTj. The only differences 
are that here we estimate the truncated value, so K has one less time integration and 
the individual integrals are performed by using Lemma I2.1L The details are left to the 
reader. □ 

Finally, if r = 1 and w G is the location of the gate/0-index among the core 

indices, then = Pw+i or = Pw +2 (depending whether we have a 0 or a gate) is 
forced by A and similarly for p^. In this case, the estimates in Lemma l5.II are worse by 
a factor of t. This factor can be absorbed into the main term . This completes the 
proof of (ITOi . □ 

5.2 Recollision and triple collision 

For the proof of we break up the partition A into the trivial partition Aq using 

Operation 1. Since s(A) < 2q, we have s(A) < 2q + 4, thus Operation I will be applied 
at most 2g + 3 times. Clearly 

|K(P(A,D))| < A2‘'+3a9sup^*(P(Ao,Do),u) (5.4) 

U 

by using and Lemma 9.5 from m- The single lumps are removed by Operation 

II from Fi*(P(Ao, Do), u), at the price A/p; the total contribution of one 0-removal is 
(AV 0 ) A The possible gates are eliminated by Operation IV at the expense of 
A^p“^| logpl ~ A”''! log A| each. Since r < I, we lose at most a factor A“^'^| logpp in this 
way (see Fig. 0 the double line denotes truncated propagators). 

The vertex set of the remaining graph is naturally identified with /c-i- 2 U/c-i- 2 U{ 0 , 0*}, 
the permutation a provides a pairing between the elements Ic +2 \ {wi,c + 2 } and Ic +2 \ 
{wi, c + 2}, furthermore {tci, c -|- 2} and {wi, c + 2} each form a lump (see picture). We 
denote this partition by P* and by using Proposition 15.21 below, we will obtain (14.351) . 
□ 

Later we need to estimate asymmtric recollision graphs as well, so we formulate the 
following proposition in a more general setup: 
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6-3k 4d 

<cx C 


Figure 7: Estimate of a two-sided recollision graph 



Figure 8: Estimate of a one-sided recollision graph 


Proposition 5.2 Consider the Feynman graph on the vertex set Vk, k > 3, choose 
numbers a,b,a',b' G Ik such that b — a > 2, b' — a' > 2. Let a be a bijection between 
Ik \ {a, b} and Ik \ {a^ b'}. Let P* be the partition on the set Ik U Ik consisting of the 
lumps {j, a{j)}, j G Ik\ {a, b} and {a, b}, {a', F} (Fig. Ql. Then 

sup E,g(P*,u) < . (5.5) 

u,3<8 

We also need a “one-sided” version of this estimate (Fig. |HJ. 

Proposition 5.3 Consider the Feynman graph on the vertex set Vk,k- 2 , k >3, choose 
numbers a,b G Ik such that b — a > 2. Let a be a bijection between Ik \ {a, b} and Ik- 2 - 
LefP* be the partition on the set IkUlk -2 consisting of the lumps {j,cr{j)}, j G Ik\{a,b} 
and{a,b}. Then 

sup .Fg(P*,u) < , (5.6) 

u,g<8 

and for the truncated version 

sup E,g{V\vL) < . (5.7) 

u,3<8 
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Proof of Provostion \5PA We use p = (pi,... ,Pk+i) and their tilde-counterparts to 
denote the momenta to express 

E^:g(P*,u)= sup [ dad/3 /3) 
g : \g\<gJ-Y 


with 


~ k 

Eg{a,l3) J d/i(p)d^(p) JJ 


\a - uj{pj) - ip\ 1/3 - u}{pj) + ip\ 

xS{pk+l - Pk+l)s(pa+l -Pa + {Pb+1 -Pb) - Ua^5(^ -Pa'+l + Pa’ “ (Pfe'+l - Pb') “ ^a') 


k 

X n ^{Pj + ^-Pj - (P’TiP + l -PcT(j)) - AAg(w) , 

37^0.,b 

where the u-momenta are labelled as u = (ui,... ,Ub-i,Ub+i,Uk,Ua’) and we used the 
identification from between the w and p, p notations. 

Without recollision, the momenta p formed a spanning set of all momenta, and 
similarly for p. Since now there is a delta function among the p momenta, we need to 
exchange one tilde-momentum (out oi Pa’,Pa’+i,Pb’,Pb’+i) with a non-tilde momentum 
(out oi pa,Pa+i,PbTPb+i)- We will call them exchange momenta. 

For the moment, we choose py and pb to be the exchange momenta and we partition 
the set of all p, p momenta into two subsets of size k + 1 each: 


A := {pi,p2,... ,Pb-i,Pb+i, ■ ■ ■Pk+i,Pb'} 


B ■■= {pi,P2, ■ ■ ■ ,Pb'-l,Pb'+l, ■ ■■Pk+i,Pb} ■ 

It is straightforward to check that all ^-momenta can be uniquely expressed in terms 
of linear combinations of the il-momenta (plus the u-momenta) and conversely. In 
particular 

Pb—l Pb {Pa{b—1)+1 Pa(b—1)') ^b 

Pb’-i = Pb' - {Pm+1 - Pm) + Um with m := a~^{b' - 1) . 

The letters on the pictures indicate the indices of the correspoding p or p momenta. 

We perform a Schwarz estimate to separate A and i3-momenta at the expense of 
squaring the propagators, but we keep the denominators with pi,pi,pb-i,pb’-iPb,Pb' 
common and only on the first power: 


k 


n 


1 

a-(D{pj) 


1 


ir]\ 1/3 - io{pj) + iri\ 


<7T[(a) + (&)] 


(5.8) 
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1 


1 


1 


(“) - n 


n 


,=1 fc-l b I" “ ^=2 1“ “ 

’ ’ j7^b-l,b 


(t):= n 


1 


1 


n 


=1 6_1 b 1“ “ ~ 3=2 

’ ’ j^b'-l,b' 


\f3-oj{pj) +ir/|2 


(with a little abuse of notations we used / for 1, b' — 1 and b' when j = 1,6 — 1 and 6, 
respectively). 

Since common factors can be explicitly expressed both in terms of A and H-momenta, 
we can compute the integral of (a) by first integrating all H-momenta that removes all 
delta functions, then estimating the ^-momentum integrals. Similar procedure works 
for (b). The result is (with m := a~^{b' — 1)) 


P p fc +1 

E*g(P*,u) < A® sup jj dadP J dp{pb')(^]^dp{pj)^Mg(w) 


3 = 1 
3#6 


\a - u;{pi) - ir]\ 


(5.9) 


1 


1 


1 


\f3 - uj{pi) + ir]\ \a - uj{pb-i) - ir?| \f3 - uj{pb' - {Pm+i - Pm) + v-m) + 

^ \2 


1 


1 


n 


A" 


\a - uj{pb+i - Pa + Pa +1 - Ua) - ipl \p - uj{pb') + irjl \a - u}{pj) - ir]\'^ 

j7^b-l,b 


We recall the key technical bound Lemma 10.8 from m to estimate integrals with 
shifted denominators. We also recall the notation |||g||| := rj + min{|g|, 1}. We start with 
estimating 

_A2_^ 

\a - w(pfe+i) - ir/|2 “la - cl;(pfc+i) - ir/| 

and integrating out p^+i by using (10.25) from [ini. We collect a point singularity 
IllPa — Pa+i + and a factor logr/p. This argument works if 6 < /c; 

the b = k case is even easier since the denominator with pk+i is not present. Then 
we perform the pb' integration again by (10.25) from ^U], and we collect a new point 
singularity |||pm+i — Pm + and a factor CC,^~^\ log Ap. Note that these two point 

singularities are not identical, because m, as an inverse image of cr, is not equal to a. 

Next we integrate out pb-i yielding C\ log A| from H2.h|) . If pb-i appears in one of the 
point singularities, then we use the bound 


sup 

Q;,r 


/ 


dh(p) 

a — uj{p) + ir]\ 


^ <CC''-2|logr?| 


(5.10) 
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(see (A. 3) and (A. 7) from [TOl) to collect log A| and the point singularity disap¬ 

pears. Ifpb-i appears in both point singularities, then we separate them by the telescopic 
estimate (A.l) of [10] before applying (Eini). 

Now we integrate out all pj’s with j ^ 1, 6 — 1, 6 , 6 -|- 1 in decreasing order with 
the successive integration scheme (10.7)-(10.9) from Section 10.1.2 of | 1 ()| . The factor 
A/’g(w) provides the necessary \B{pj — pj_i)p terms with at most eight exceptions, 
namely when j — 1 G Q oi a{j — 1) G Q (recall |^| < 8 ). At each exceptional index 
\B{pj — is replaced with {pj — and we use (j2.7l) instead of (10.8) of 

m- Thus the successive production of the factors (1 -|- C\^ breaks at these indices 
and we obtain a uniform constant C instead. The successive scheme also breaks at the 
indices j = b — l,b,b + 1 that have already been integrated out. Furthermore, it also 
may break at j = m -|- l,a -|- 1, i.e., at the indices where the point singularities are 
first affected (unless 6 — 1 G {m -|- 1, a -|- 1} and the point singularity has already been 
integrated out). At each of these indices we use H5.10|) and collect CX‘^r]~^C'^~‘^\logri\‘^ 
instead of the constant factor from (10.7)-(10.9) of |10j . 

Since there are at most 13 exceptional indices, so we collect at most 

logr/p]2(l + CX^-^^'^)^. 

The other factors of Mg, that are not explicitly used in the successive integration, are 
estimated by supremum norm, except |V'o( 7 'i)P- Finally the da, d/9 integrals contribute 
with an additional CIlogAI^. The last pi-integral is finite by the factor Col¬ 

lecting these estimates and recalling that (lOH) has been used twice, the result is (ESI). 

□ 


Proof of Proposition \5JA This proof is very similar to the previous one but the 
estimate is weaker since (10.25) from (TUI can be used only once. We just indicate that 
the set of A and B momenta are as follows: 

A := {pi,P2, ■ ■ ■ ,Pb-i,Pb+i, • • • ,Pk+i] , B := {pi,p2, ■ ■ ■ ,Pk-i,Pb} , 
and we leave the details to the reader. □ 

The case of the triple collision, can be easily reduced to the case of a recollision 

(Fig. El). We first use the analogue of the estimate dOl). Clearly <7 = 0 in case of a triple 
collision. Then we remove half of each of the two gates (Lemma EiSl) and we collect a 
factor A^l log Ap. The resulting Feynman graph has either a recollision or a gate at the 
end. In the first case we apply Proposition 15.21 In the second case, we remove half of 
each gates by a second application of ljemma. lC.3l then the estimate ()C.2|1 from Lemma 
ICJ together with (EH) can be applied. □ 
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Figure 9: Estimate of triple collisions for w < c+1 and tc = c + 1 

5.3 Cancellation with a gate 

The key mechanism behind the estimates is the cancellation between a 

gate and a 0-label, We first present estimates on general graphs. 


5.3.1 Cancellation between a gate and 9 

Fix n, n' integers and consider a partition P £ Vn,n' with no single lump on the set 
In U In' within the vertex set V = Vn,n’- Let l<m<n-|-lan integer. We define two 
new cyclically ordered sets: 

VC= {0,1,2,... ,m — 1, ... ,n,0*,h',n' — 1,..., 1} (5-11) 

V := { 0 ,1, 2 ,... ,m - 1, 0 ,^,m,... ,n, 0 *,h',n' - 1 ,..., 1} 

with additional elements Jlk, <C>, For the result of this section it would make no differ¬ 
ence if the extra elements were inserted into the sequence of tilde-variables. 

These sets can be naturally identified with Vn-i-i,n' and Vn+ 2 ,n’ and we will use this 
identification with the obvious choice of the relabelling map. We define two partitions 
on these sets, P' £ 'Pn-i-i,n' and P" £ T’n+ 2 ,n', simply by adding the single lump {JH} to 
P in the first case and the double lump {0, *7} in the second case. This will correspond 
to adding a 9 label or a gate whose potential labels have been paired to the original 
partition P, respectively. The following lemma shows that the F-value of these two 
partitions cancel each other up to the lowest order fFig. 11011. 


Lemma 5.4 With the notations above and assuming rj we have 


F(,)(P') + F(,)(P") 




(5.12) 
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Figure 10: Cancellation of a gate and 6 (partition is the same elsewhere) 


Proof of Lemma\^^ Introduce the notations p := (pi,... ,p„_|_i), p := (pi,... ,Pn'+i), 

j dp := j j dpidp2 ... dpn+i 

and similarly for f dp. Then we have 

p2tr/ p pY p 

Fp)(P')+F(,)(P") = j J dadd y dp dp A(P, w, u = O)M(w) 


n+l, (n) 


n'+l, (n') 


:Q{a,pm) 


f=i oi - uj{pj) - ir] P-uj{pj)+ip' 


n 


with 




\B{q - Pm^dq - 


a — uj{q) — ip 


: -6'(Pm) 


A 2 


a - U){pm) - ip 


(5.13) 


(5.14) 


The expression n + l, (n) on the product sign indicates that for the truncated values 
(*) the last fraction is not present, i.e. j runs up to n. Notice that the sum of the 
contributions of a gate and a inserted between m and m — 1 yields an additional factor 
fl(a,Pm) in the F-value of the P partition. Using (|B.1|) and (|B.9I) with e = 77 , e' ^ 0 + 0, 
we have 

[ _ e{pm) < + p~^/‘^\a - A^e(a) - e(pm)l) • 

J a — uj[qj — ip \ / 

Therefore, using (IB.5j) and <C p <C A^, we have 

|S!(a,P„)l < CAV'hl + 1!!++:^) < (5.15) 

V \a-UJ[Pm)+ip\ ^ 

uniformly in a and Pm- d 


The proof shows that the cancellation between the gate and the i? is completely local 
in the graph. In particular, if we fix L locations, maybe with multiplicity, between the 
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elements of and we consider all possible 2^ combinations of insertions of gates and 
7?’s at these locations, then we gain a factor from each location. 

More precisely, let v € be a given sequence of integers, vq,vi, ..., vq* , labelled 

by the elements of Vn,n'- The number Vj indicates how many gates or i?’s are inserted 
between the j-th and (j —l)-th vertex. Let |u| := Vj be the total number of insertions. 
A sequence S G encodes whether the insertion is gate or -d. 

Fix a sequence v and for any S G {g, we dehne the extended set V 5 consisting of 
Vn,n' and we insert extra single or double symbols for i} and gate indices (determined by 
S) at the locations given by v. In the example above, we have v = (0,0,..., 0,1,0 ... 0) 
(i.e. = 1, the rest is zero), |u| = 1, and V' corresponds to 5 = {?9}, while V" 

corresponds to S' = {^f}. Given a partition P G Pn.nG we also dehne the extended 
partition P5 on V5 by simply adding the single symbols as single lumps and the double 
symbols as paired lumps to P. 


Lemma 5.5 With the notations above, for any fixed P G Vn,n' o.nd v G 




<(caV'/')'"'^w(p)- 


(5.16) 


Proof. Notice that the sum of the contributions of a gate and a 'd inserted at the same 
place between m and m — 1 yields a factor of Q{a,pm), while the same insertion between 
fh and m — 1 yields a factor Q(/3,pm). These insertions are independent of each other, 
thanks to the summation over all possible 5-combinations. So F(P 5 ) is represented 
by an expression similar to Km . where a total factor 

[n(a,Pm)]‘'"‘ [n{P,PniW^ 

m=n' ,...,1,0 

is inserted. The uniform estimate for each factor gives (tPriHl) . □. 

Next we will apply Lemmas I5.4H5.5I to prove (l4.,‘-{7ll - (l4.59jl . The difficulty is that 
these estimates hold only if the gate remains isolated even after the lumping procedure, 
otherwise the gain comes from the artificial recollision introduced by the lump. Recall 
that the lumping procedure has two steps. The original partition Dq may lump non- 
trivially, yielding the derived partition D, due to a few possible coincidences among the 
gate or recollision labels of and fi. Then the non-single elements of D (denoted by 
D*) lump into a coarser partition imposed by A G A.(D*) due to the connected graph 
formula. 

We will estimate the value of individual graphs only. The number of terms in the 
summations in (I4.d7l) - (l4.;-i9l) is bounded by 0((c -|- 4)^) < CK'^. This extra factor 
will be added to the factors gained in the cases discussed below to obtain (I07ll - (I0^ . 
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5.3.2 Non-repetition graphs with a gate 

To prove (1071) . we first note that the non-repetition rules in force D to be 

identical with Dq unless h, h' = g. In this latter case there is a gate both in the expansion 
of iIj and and either D = Dq or D lumps the two gate-lumps in Dq together. 

We first consider the case D ^ Dq. Then D has a lump consisting of all four gate 
indices, {re, w + 1 , 10,10 + 1}. The corresponding Feynman graphs can be identified with 
certain nontrivial lumpings of non-repetition graphs on Ic+2 U/c+2, where the indices w 
and w + 1 are lumped together. More precisely, for a given a e ©c, w ^ Ic, h = h' = g, 
G Ac, 

F(P(A,D))c(A)= ^ V{A,a)c{A), (5.17) 

Dt^Dq AeA.(D*) Ae-4,c-i-2vA— 

w=w+l{mod A) 

where a G ©c-1-2 is the natural extension of a where two new elements, w,w + 1, are 
added to the base set Ic, the indices are shifted by the embedding map s!^ fSection l4.3|) 
and a{w) = w, a{w -|- 1) = w + 1. The summation has at most c -I- 1 terms and it 
expresses the choice of joining w,w + 1 to one of the existing lumps in A' or keeping 
the lump {w,w + 1} separate in A. Because of lumping w and w + 1, the partition A 
is non-trivial, and q[A, a) > 1 (see ). After estimating |F(- • • )| by E{- ■ ■), each 

term on the right hand is estimated by the bound 

suvE^,)g{A,a,u) (5.18) 

U ^ ^ 

that holds whenever cr G ©fc, A G Ak and n < 34 ^^ (see (9.4) from [Ej)- 

The estimate (EHI), combined with the combinatorial bound on c(A) and on the 
summation (Lemma ESJ; gives (I4..‘17jl in the case when D A Dq- 

Now we focus on the case D = Dq. If at least one of the gate-lumps, {w,w + 1} or 
{w,w -|- 1}, do not remain isolated in A, then we can repeat the argument above since 
P(A,D) has a non-trivial lump of size at least 4. 

Finally, we can assume that the gate lumps remain isolated in A. We fix A' and 
consider the sum of four terms corresponding to h, h' G {g, 0}. The partition A is defined 
by adding the gate lump(s) to A'. Note that c(A) is the same for all these four cases 
since replacing a 9 index with an isolated gate lump adds only a single lump to A. For 
these partitions, we apply Lemma [531 fwith the choice Vw = v^, all other u’s are zero) to 
obtain a factor (CA^“ 2 )^. The remaining non-repetition graphs bounded by 0{\ log Ap) 
by using Proposition 9.2 from m with q = 0. This completes the proof of (lOTll . □ 
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Figure 11: Case 1. Removal of a gate lumped to an adjacent core index 

5.3.3 Last gate 

First we consider the case when the lumps {re, tc + 1} C In and {w, tc + 1} C In' of the 
two hrst gates remain isolated in P = P(A,D). As in the previous section, by applying 
LemmaESl we can sum up the two times two possibilities for the first components of the 
code h and h', i.e. sum up four Feynman diagrams that differ only by the choice of gate or 
9 at the w ox w position. We collect The remaining two gates will be isolated 

from the rest in P by Operation I, then half of each gate is removed by Operation III, 
collecting A^A^| log Ap. Finally, (I0.2|l can be used to remove the remaining two halves 
of the gates, collecting A^. By (15.Ij) . the untruncated Fi-values of the remaining graph 
are bounded by 0(| log Ap). The total estimate is 

If only one of the lumps, {w, tc +1} or {ui, re + 1}, remains isolated in P, we can still 
apply Lemma El to obtain a cancellation of order A^ from adding up the F-values 
of those pairs of graphs that differ only by changing this gate to 9. Note that the value 
c(A) is again the same for these two graph. 

Next we consider those lumps among {w, tc + 1} and {w, rc + 1} that do not remain 
isolated in P. By breaking up lumps via Operation I, we can ensure that every such gate 
is either lumped exactly with one other gate or with a core index pair. For definiteness, 
let {re, rc +1} from be such a gate. We also isolate all other gates from the rest. To do 
that. Operation I is used at most four times at the total expense of A^. We distinguish 
three cases: 

Case 1. The gate {tc, rcT 1} is lumped with a core index-pair (j, a{j)), where a is the 
natural extension of a from Ic to In- If j is next to w or w+1, say j = w-\-2 fFig. lllI) . then 
both vertices of the gate can be removed by Operation III since the momenta between 
{w — l,tc) and {'w,w -|- 1) do not appear in any delta function, and we gain A^| log??p 
from this gate. We can now remove the two gates at the end (using Operation III and 
as above), collect an additional A"^! log Ap. The remaining gate at {w,w + 1} can 
be removed by Operation IV at the expense of A“'^| log A|. By Proposition 9.2 from | 1 ()) . 
the remaining graphs are bounded by 0{\ logAp). We thus collect C 

If j is not next to u) or tc +1 ("Fig. 11211 . then we remove tc +1 by Operation III, break 
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Figure 12: Case 1. Removal of a gate lumped to a non-adjacent core index 



Figure 13: Case 2. Removal of two adjacent lumped gates 

up the lump into {j, rc} and {?(j)} by Operation I and remove the single lump 
by Operation II. The total price for these steps is logr/p. We can now again 

remove the two gates at the end, collect A'^l logAp and we end up with a graph with a 
one sided recollision, so ((SH) from Proposition Eini applies. We collect CA® (i 4 rf+ 2 )K o(<5) 
in this case. 

Case 2. The gate {tc, u) + 1} is lumped with the other gate in I„, i.e. with {n — 1, n}. 
If u) + 1 and n — 1 are neighbors (Fig. 1131) , then we can remove three vertices n — 2 = 
tc + 1, n — 1, re by Operation III. We collect A^| logr/p. On the In' side, we remove the 
gate that is not adjacent to 0* by Operation IV at the expense of A”''! log A| as above. 
The other gate is adjacent to 0*; first we remove its leg not adjacent to 0* by Operation 
III, collecting A|logA|. Finally, we remove the two remaining vertices, originally with 
indices re — 3 and n' that now form single lumps and they are both adjacent to 0*. By 
using (EI21), we gain a factor A^. Altogether we thus collect CA® 0 ( 5 )^ 

If rc + 1 and re — 1 are not neighbors iFig. I14|l . then we remove tc + 1 and re, gaining 
A^|log? 7 p and the remaining partition again has a one-sided recollision. The other 
gate adjacent to 0 * can be removed by collecting A^| log Ap, the remaining gate collects 
A“'^| log A|, and finally we use (|3.(ij) . The result is again ^ 6 -(i 2 (i+i)K- 0 ( 5 ) _ 

Case 3. Finally, if neither {tc, tc -|- 1} nor {rei, re -|- 1} falls into Case 1 or 2, then each 
of them either remains isolated and collects from Lemma E31 or is lumped with 

another gate on the “opposite side” (Fig. 1131) . For definiteness, assume {w^w + 1} is 
lumped with 1} in /„/, where (. is either tc or re' — 1. Then we simply remove half 
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Figure 14: Case 2. Removal of half of each gates lumped on the same side 
g g+i 



Figure 15: Case 3. Removal of two opposite lumped gates 

of each gates, say re + 1 and £ + 1 using Operation III, gain A^| logryp and we extend 
the set of core indices to include w and extend the permutation a by adding a{w) = i. 
Therefore we effectively gained A| logr;| from each such gate. Finally, after having gained 
either or A| log A| from each gate, we gain A^| log Ap from the remaining truncated 

graph (I5.1|) -(I5.21) and we thus collect at least CA® 0{5) ^ This completes the proof 

of (lOHl) . □. 

5.3.4 Nest 

The procedure is very similar to the analysis of the last gate, so we just outline the steps. 
We point out that the main reason why the nested graphs are small is the cancellation 
between the gate and 6 inside the nest. This is a different mechanism than the one used 

in [ 3 . 

First we consider the cases when both nests are independent, i.e. no part of the 
nest in is lumped with any part of the nest in /„/. We will show that one can gain 

at least ^3-(2d+2)K-0(<5) each such nest. The total gain from both nests is then 

^6—(4(i+4)K—0(5) 

Consider the gate {n — 2, n — 1} inside the nest. If this gate remains an isolated 
lump in P, then it cancels the same graph with 6 up to order by Lemma 15.41 

After this cancellation, the outer shell of the nest, {n — 3,n}, becomes a gate in the 
reduced graph that is adjacent with 0*. After separating it from the rest by Operation I, 
if necessary, its removal yields A^ because of the truncation. Therefore we gain at least 
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AA3-K/2 < ^3-(2d+i)K-0(5) 

If the gate {n — 2, n — 1} is not isolated in P, then we again distinguish a two cases. 

If {n — 2, n — 1} is lumped with a core index j < n — 3 but not with its outer shell, 
{n — 3, n}, then we can always create a one-sided recollision in such a way that n — 2 will 
be paired with a{j) in the extended permutation, while n — 1 is removed (Operation III, 
gain A| log A|) and j is removed (Operation II, lose The net result is A“'^| logAj 

and the outer shell of the nest, {n — 3, n}, becomes a genuine recollision. We may have 
to separate this from the rest of the graph by an Operation I before applying a. 
This will effectively give AA^ log Ap^^^ < A^ { 6 d+ 2 )K o{5) £qj. 

calculation we used only A^”'' from (|A.7II . the additional A is due to the truncation on 
the “other side” and will be counted there. 

If the gate {n — 2, n — 1} is lumped with its outer shell, then the momenta between 
(n —3,n — 2 ), (n — 2 ,n —1) and (n —l,n) can be freely integrated, we can remove the nest 
completely and collect AA'^ < CA‘^“^'^'^ from this nest. We may have to use Operation I 
once to separate the nest from the rest of the graph. 

So far we have treated the cases when the two nests are independent. In the remaining 
cases some parts of the nests are lumped with each other. 

If the gate {n — 2, re — 1} is lumped with the other gate {re — 2, re — 1}, then half 
of each gate is removed, say re — l,re — 1 by Operation III (gain A^| logAp), the other 
halves are separated from the rest of the graph (Operation I) and then connected by 
an extending the permutation cr(re — 2 ) = re — 2 (i.e. we include re — 2 among the core 
indices). The resulting graph has two recollisions, so (15.51) applies and the total size is 
AA®-3«C4'^|logAp(i). 

Finally, if {re — 2, re — 1} is lumped with the outer shell {re — 3, re} of the other 
nest, then we isolate and remove the gate {re — 2 , re — 1 } inside the other nest (price: 
AA”'^! log A|), remove re and re — 1 each as one half of a gate, gaining A^ and the remaining 
graph is a one-sided truncated recollision graph after a possibly isolating the recollision 
lump {re — 3, re} from the rest. Thus (EH) gives Aa 4 jQg ^| 0 (i)^ after a possible 
application of Operation I. The total size is Collecting the various 

cases, we obtain (ICTll . □ 

6 Convergence of the ladder diagrams to the heat 
equation 

We start the proof of Theorem 12.41 by noticing that 

wx{t,k,o) = J v4(Ao,a(e,-))de, k>i 
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with Aq being the trivial partition on where we chose the function Q{v) in the 
definition of V° to be ^-dependent, namely Q{v) = Q^{v) := 0{^,v) (see (j4.17l) for the 
definition of V°). First we note that the d^ integration can be restricted to the regime 
{I'^l < with a negligible error (even after summation over k): 

Wx{t,k,0)= + /'V;^(Ao,a(e,-))d^, (6.1) 

l<k<K l<k<K 

where use the notation 

f {■■■H ■■= I {■■■)m\<x-')dc. 

To see (EO, we first recall that replacing 1^°(- • •) with !/(• • •) yields a negligible error 
even after summing up for all k (Lemma 7.1 of ^HI). The linearity of the estimate in 
IIQ^Iloo = sup^ |C>(^,u)| guarantees the integrability in ^ since O is a Schwarz function. 
We then use the estimate 


14 ^(Ao,a(|,-)) < ||Q^||ooSupF;(cr = id,u) 




and the uniform bound m and finally we conclude (iniii) by the arbitrarily fast decay 
of IIQgIloo in 6 

From the definition of we have 


II* de I n Oi^, Vk+i)WoH, ^i) 

k-\-l "p ^ 

X Yl - Y)\B{vj - Vj+i)\'^ . 

i=i 


( 6 . 2 ) 


To simplify the notation, in (jnn we followed the convention, that \B{vj — Uj+i)| = 1 for 
j = k + 1 because of the non-existence of Vk+ 2 - Similar convention will be followed later, 
also for \B{vj-i — Vj)\‘^ = 1 if j = 1. Note also that the measure dvj has been changed 
to dfj,{vj) by using the support properties of B and ipo (12.11) . 

The estimates of the error terms were performed with the choice rj = However, 

given by (lO) . is clearly independent of rj; this follows from the A-identity (formula 
(6.2) in [1 Dji. Therefore we can change the value of ry to r/ ;= for the rest of this 

calculation and we define 


R{a,v) := Rri{a,v) , with ry := . 
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We also recall, that the restriction of the dad/9 integration in (inii to any set that 
contains {a,/9 : |a|, |/9| < y = results in negligible errors, even after the sum¬ 

mation over k (Lemma 7.1 of [I^). We will consider the set D := {(a,/9) : \a +13\ < 
2Y, \ a — f5\ < 2Y }. We denote by the version of given by formula (16.21) with the 
dad/9 integrals restricted to D, 




\2k f f dad/9 

JJnWr 


Integrand from (ESI) , 


then 

l<k<K 

We also remind the reader that this argument literally does not apply to the trivial 
k = 0 case, when the da d/9 integral in (ESI) gives free evolutions and this term should 
be computed directly: 

“0 := j d^du d{^,v)Wo{e^,v) + o(l) , (6.3) 


where the error term comes from the error term in 6 (v + £^/ 2 ) — 6 {v — £^/2) = 2 il{y) + 
0{e^). By using oo, the bound 

Im 9{v) < —Cl min{|p|'^“^, |p|~^} , 


(from Lemma 3.2 of [TO] ! and the decay of the observable, one easily obtains that |Ho| = 
o(l) anyway. 

To evaluate the integral EH, we need the following crucial technical lemma which 
is proven in the Appendix. 


Lemma 6.1 Let k < 1/8, define 7 := (a + /9)/2 and let rj satisfy < r] < 

Then for |r| < we have, 


n : = 


AV(p) 


^a — u}{p — r) — irj^ ^/9 — uj{p + r) + 


dp 


= -2-KiY 


f{p) 6{e{p) - 7 ) 


{a — P) + 2p ■ r — 2i[\‘^Z{p) + p] 


dp + 0{\^/‘^ 


-4k 


4d,l 


(6.4) 


Now we compute by applying Lemma I6.1L Denote a := (a -|- /9)/2 and b := 
A“^(a — /9). By using e = p = we have 


A^ / dvT{f,v) R{a,v + 


V — 




(6.5) 


-2mT{i,v) 5{e{v) - a) ^ _ i/ 2 - 4 kmi^ii 
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In the applications, T(^,u) will always be supported on |ri| < (, therefore the measure 
dr; can be freely changed to d^(r;). 

We now replace the product of /c + 1 factors in the restricted version of one 
by one. We need a factor for each application of thus we need in H(i.2l) . 

But contains only the missing A^ comes from the change of variables dad/3 = 
A^dadfe. The domain of integration, (a, /3) G D, is replaced by the domain D* := {(a, b) : 
|a| < F,|6| < 2A-2y}. 

For any £ = 1,2,... k + 1, we introduce the notation 


k+l 


£-2 




J =1 


27ri Fj(^,vj)lB(vj S(e(vj) - a) 


b + ^ — 2i[T(o) + A^'^] 


2TriFi{^,Vi-i) 5{e{vi-i) - 

-a) 

b + X^/^Vi-i ■ C - 2i[T(a) + A^^] 

j diJ.{vi)Ti{vi_i,Vi,Vi+i) 

1- 

> 

to 

^ 1 


X^R(a,Vi + - y) 




—27rz 6{e{vi) — a) 


b + C - 2i[I{a) + 


X n [>f 

j=e+i \ 






R[a,Vj + —)R[f3,Vj -^)Fj{^,Vj) \B{vj - Vj+i)\ 


X A^ 


R(^a,Vk+i + y)i?(/3,Ufc+i - y^Ffc+i(^,Ufc+i) 




( 6 . 6 ) 


with 


Ti{ve_i,Vi,Vi+i) 


Wo{e^,vi) U 2 )P for £ = 1 

< - Vi)\^\B{vi - Vi+i)\‘^ for 2<£<k 

|0(^,r>fc+i)| |.B(ui - U 2 )P for £ = k + l 


and 




VFo(e^,u) for j = l 
< 1 lov 2 < j < k 

0{^, v) for j = k + l ■ 


The formula ( 110 ) 1 ) is literally valid for 2 < £ < k. For £ = 1 the first product and 
the factor in the second line are absent, for £ = k + 1 the factors in the last two lines 
are absent. We also recall the convention made after (lO) about the intepretation of 
\B{vj - Uj+i)| for j = k + l. 
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With these notations and by introducing r := = A and W{k) := W\{t, k, O), 

we obtain the following telescopic estimate from (lini) 


^W(k) 

k<K 


/ * 

d{ 


'D* 


dadb irb+2tri 
( 2 vr )2 


X 


( ti [ -27rzF(^')(.eWi) Kejvj) - «) I 
V J =1 J b + X^/^Vj ■ C - 2i[I{a) + A^'^]' ^ ^ 


fc+i 

< ^^k,£ + 0 ( 1 ) . 

k<K £=l 


Vj+i)\‘^dn{Vj] 


(6.7) 


Now we explain how to estimate J-k/ for the general case {2 < i < k), the modifications 
for the two extrema are straighforward. 

First we estimate Wq by supremum norm and estimate all denominators in the first 
two lines by their imaginary part: 


1 

b + ^ - 2i[Z{a) + A4«] 


1 

“ 21{a) ■ 


( 6 . 8 ) 


Then the ui, U 2 ,..., U £_2 variables are integrated out in this order, by using JH, yielding 
a total factor 1 from the product in the first line of (IHIH) . By recalling (tOl) and the 
estimate 

X(a) = —Im 0(a) > ci min{|e| 5“^, 
from Lemma 3.2 of m, the integral of U£_i is estimated trivially by 


1 

21{a) 



a) di 2 {vi-i) < 


Ca^/^ 

I{a) 


< (a) . 


(6.9) 


This estimate is used if a < C^/2, otherwise the integral is zero by the support of d/r, so we 
obtain a factor In the regime \b\ > A“'^, we have \b — • .^| > \b\/2 

using \vi-i\ < C and |^| < A””^. The estimate (j6.8l) can thus be changed to 2\b\~^, 
improving estimate dOH) to < CaV2/|6| < CC/\h\. 

The integral d^{vii) in (16.61) is estimated by 

0(A1/2-4-) sup ||T^(u^-l, -,^^£+1)114^1 < 0(A^/2-4K^4d) ^ Q(^l/2-(4d+4)K-0(5)) 

( 6 . 10 ) 

by using (16.5() and the fact that all vj variables satisfy \vj\ < C- For ^ = 1 and £ = A: + 1 
we also used that the initial data and the observable are Schwarz functions. 
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For the dfj,{vj), j = £+l,i + 2,... ,k, integrals we separate the resolvents by Schwarz 
inequality, 

k k k 

\R{a,Vj + . ■ ■)R{(5,Vj - .. .)\ < \R{a,Vj + .. .)\‘^ + \R{P,Vj - .. , 

j=e+i j=£+i 3=^+1 

and we use the successive integration scheme (see Section 10.1.2 of cni) to collect a 
constant factor. 

Before we integrate out the last momentum variable, Vk+i-, we perform the dad 6 
integration. We can change back the a,b variables to a, (3, we perform dad/3 integrals 
to collect a C|logAp factor since D C {|a|,|/3| < 2Y}. This argument applies unless 
i = k + 1 and the last line in is absent. In this case, however, i > 2 (the k = 0 
case is treated separately, see (EH), and then the denominator with U£_i in the second 
line of (10.61) is present. We use the a, b variables. Recall that 5{e{v£-i) — a) restricts the 
domain of the da integration to \a\ < C^/2, giving a contribution O(C^). The domain of 
the b integral is larger, | 6 | < 2A“^y, but in the regime |5| > A“'^ we have collected an 
additional factor in (16.9j) . thus the db integration contributes at most with a factor 
0 (A-^). 

Finally we integrate Vk+i by using the integrability of T^+i = O and the ^ integral 
gives a factor 0(A“^^). By collecting these estimates, we arrive at 

fc + l 

E E 

l<k<K £=1 

and that is negligible, since k < l/{8d + 18). 


Now we focus on the main term on the left hand side of dnizi)- First we extend the 
db integration from |6| < 2A“^y to M. It is easy to see that the error is negligible; all 
denominators can be bounded by \b\/2 and the result from the db integral in the region 
|6| > 2A-2y, 

f < X2k 

Jm>2X-^Y - ’ 

is negligible even after multiplying the from the dvj integrals. We can also extend 
the da integration from [—Y, Y] to M, since, due to the factor S{e{vj) — a) and the cutoff 
in Vj, the integrand is zero for |o| >Y. 

Now we write 


—i 

b + C - 2i[I{a) + A4«] 
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and perform the db integration. We obtain 


/ * r k-\-l pQQ fc+i 

^hXr-^r,)e 

j = l j = l 


-iA«/2(52 


fc+l 


^ (n / 27r|-B(z;j - z;j+i)p5(e('(;j) - a)j6(^,?;A:+i)TVo(eC,^^i) + o(l) • 

i=i 


We now replace d/x(t'j) with and remove the cutoff in ^ to perform the Fourier trans¬ 
form. We also replace Wq with Fq and remove the k < K cutoff from the summation: 


Y,wik) 

k<K 


CXD „ „ k + 1 „oo ^ + 1 

k=0'^ j=l 


/ k-\-l 


X J dui 6{e{vi) — a) ^ JJ J dvj 2Tr\B{vj — Vj-i)\‘^6{e{vj) 
X 0{X,Vk+i)Fo(^X - {27r)-^X^/^C^TjVj), ui ) + o(l) 


a) 


( 6 . 11 ) 


with initial data 

Fo{X,v) :=6{X)\Mv)f ■ 

The replacement of dfj,{vj) with dvj is justified since V’o('i’i) is compactly supported 
and thus all other Vj’s are restricted to a compact energy range by the delta functions 
Wj5{e{vj) — a). The removal of the ^-cutoff is allowed since the integrand of the d.^- 
integral can be majorized by 


sup|0(^,u)| f 2 tT(a) ^ [2X(a)f] f s{e{vi) - a)Wo{e(,vi) < sup\d{^,v)\ 

V Jr 27r ^ k\ J v 

^ ( 6 . 12 ) 

whose integral vanishes in the regime |^| > A ^ due to the assumptions on O. Here we 
used (E31) to perform the Vj integrations successively and the time integration yielded 

T^/k\. The replacement of VFo(e^,u) with \'tpo{v)\'^ comes from the uniformly integrable 
bound and from the uniformity of the limit ||V'o(' £?) ~ V’o(')ll —> 0 as ^ runs 

over any compact set. Finally, the removal of the k < K cutoff in the summation follows 
from the same majorization as together with 


E 

k>K 


[21{a)r]^ 

k\ 


< 


jCr)^ 

K\ 


< {CX^)^ 0 . 
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For a fixed energy e > 0 we consider the continuous time Markov process {n(t)}t>o 
on the energy surface Eg with generator This process is exponentially mixing 

with the uniform measure on Eg being the unique invariant measure (see Lemma lA.il in 
the Appendix). Let Sf denote the expectation value with respect to this process starting 
from the initial state ip = ipo given by the normalized measure 

\'ip{v)\‘^ 5(e(ni) — e) dn 

PFIL) 

on Eg (for the notations, see Section P). Let Tg denote the expectation with respect to 
the equilibrium. Let dj 2 ^{e) = [|'!/;p] (e)de be the energy distribution of pj. The coarea 
formula, 

/ CO ^ r ^ 

[|^P](e)de = J \i>{v)fdv 

and Ip £ L'^ guarantee that dfj,^ is absolutely continuous. 

From and r = A we have 

poo pr 1 

''^^W{k)= / T^C>(A'^/^x(r), ^(t)) d/r^(e) + o(l) with x(r) := / —n(s)ds . 

k<K Jo Jo 27r 

Due to the exponential mixing and the continuity of O, the replacement of with the 
equilibrium measure Tg gives a negligible error since r ^ oo. By the central limit theorem 
for additive functionals of exponentially mixing Markov chains, converges to a 

centered Gaussian random variable with covariance matrix 

Tg[A^^^x(r) (g) A^^^x(r)] = ^2 [[ Tg [u(s) (g) n(s')] dsds'—> 2rD(e) . 

(27r)"= J Jq 

Since the equilibrium measure is uniform, the covariance matrix is diagonal, Dij{e) = 
De6ij. The diffusion coefficient, Dg, is hnite and positive. This proves Theorem 12.41 □ 

A Mixing properties of the Boltzmann generator 

The Boltzmann velocity process with generator Lg introduced in (fTT^ enjoys very good 
statistical properties. The proof uses standard arguments which we only indicate below. 

Lemma A.l For each e > 0 the Markov process {u(t)}t>o with generator Lg is uni¬ 
formly exponentially mixing. The unique invariant measure is the uniform distribution, 
[• ](e)/[l](e), on the energy surface Eg. 
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Sketch of the proof. Let V^{u,A) be the transition kernel for any u G Sg, A C Sg. 
Since the transition rate a{u, v) is continuous on Sg and 0 G supp (B) holds, there exists 
an open set S C diam(S') < and there exists and 6 = 6{e) > 0 such that 
a{u, v) > 5 whenever u — v £ S, u,v G Eg. Since the state space Eg is compact it follows 
that the transition kernel satisfies a uniform Doeblin-type condition, 

inf [ V\u,A) dt>C{e)\A\ kl C Eg , 
neSe Jo 

with some e-dependent positive constant, where |A| is the restriction of the Lebesgue 
measure (on Eg) of the set A. It is clear that the Markov process {'e(t)}t>o is irreducible 
and aperiodic, therefore it is uniformly exponentially mixing. Moreover, the rate of the 
mixing is uniform as e runs through a compact energy interval since in this case C{e) is 
uniformly separated away from zero. It is easy to see that the uniform measure on Eg 
is invariant and by exponential mixing it is the only invariant measure. □ 


B Estimates on Propagators 

B.l Proof of Lemma 12.IL 

The following lemma proves (EH) and (|2.7I) . The proof of (12.8f) is analogous but easier 
and will be omitted. 


Lemma B.l Let n < 1/6 and rj satisfying < f] < For any 0 < a < 1 we 

have the following approximation result 


a — u){p) + ig a — e{p) — A20(q;) + ig 

\h{p-q)\ 


2—a 


< c\ 




\h{p-q)\dp 
dp . 


\a — e{p) — A20(a) -|- ig\‘^ 

Moreover, for any 0 < a < 1 we have 

^ CJ|h||2d,oA-2(l-“) 


\h{p-q)\ dp 


and 


l« - e(p) - A20(a)+fr/|2 “ {\q\ - y/^\) 

\h{p-q)\ dp C'all/ilbd.o |logA| log(a) 


a - e{p) - \^Q{a) + ig\ {a)^l‘^{\q\ - ^/^\) 


(B.l) 


(B.2) 


(B.3) 
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Proof. To prove dnu) , we rewrite it as 

1 


1 


a — u){p) +17] a — e{p) — A20(a) + ir] 
A2(0(e(p))-0(a)) 


2—a 


\h{p-q)\ dp 


2—a 


(B.4) 


(a — uj{p) + irj){a — e{p) — A^0(q:) + ip) 
From the Holder continuity, 

|0(a)-0(a')l <C\a-a'\^/‘^ 
(Lemma 3.1 from |lU|i. we can bound (IH.4I) by 

|e(p)—A^ 2-a 


\h{p - g)| dp . 


(B.5) 


(a — oj{p) + ip) {a — e{p) — A20(a) + ip) 


\h{p - q)\dp . 


Since 


\e{p) — < \oj{p) — + 0(A) , 

this integral is bounded by 

\h{p-q)\ dp 


{\^p ^/‘^+\^p 


a — e{p) — \^Q{a) + ip\^~ 


Notice that in these estimates we used that the imaginary part of U7{p) is negative. 

To prove the estimates (IB.2I1 and (lESl), we rewrite the integrals by the co-area 
formula (0 < a < 1): 


\h{p-q)\dp 


f 


ds 




with 


\a — e{p) — X‘^Q{a) + ip\^ “ Jo \oi — s — X^Q^a) + ipl"^ “ 
\h{p-q)\dp 1 f 


(B.6) 


S(s) : = 


'\p\=V^ |Ve(p)| J\p\=VTs 

Using the decay properties of h, we have 

\2d,0 y/s 


\h{p-q)\dp 


E{s)< 


- VTs) is) ’ 


so 


\h{p-q)\dp 


\a — e{p) — A20(a) -|- 


< 


C\\hhd,o r 
Jo 


^/s ds 


(•s) (kl “ V^)\ot — s — X‘^Q{a) + i? 7 P “ 
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For a = 1 the last integral can be directly estimated by C | log r/] (a)' log(a), 

yielding (ESI). 

To prove (El, i-e. for a < 1, we recall that that Q = lZ — iI with non-negative real 
X and TZ. The last integral is estimated by 


f 


y/s ds 


0 (s) {\q\ - v^) [{a - s)2 + (A2J(5) + r?)2] 

with a := a — X^TZ{a). We used the Holder continuity of X, \^X{a) = A^T(5) + O(A^) 
and the fact that the error can be absorbed into rj. 

First we assume that \a\ < 1, then the estimate on 0(a) from Lemma 3.2. of |1()] 
yields X(a) < ci|a|^/2, xhe I integral can be estimated 


I < 


y/s ds 


{\q\) Jq [(5_s)2 + (A2|5|1/2+^)2]W2 +72 {s) {\q\ - V^s)S^-'^ 


i: 


y/s ds 


The second term is bounded by (|g|) ^ (l9l - a|) In the first term we consider 

two cases. If |a| < /? := A^|a|^/^ + r], then 


I 


y/s ds 


0 [(a- s)2 + (A2|a|V2 +^) 2]1 “/2 ^J\a-s\<2P J\5-s\>2f3 


</3-(2-“) / _^ds + 
J\s\<2/3 


ds 


|a—s|>/3 1*^ 


-s|3/2 


_ = 0(/?“-V2) < Q(^-2(l-a)^ 


by using that in the second regime s and |s — a| are comparable and that /3 > r/ > A^ i 
the last step. This proves El for |a| < 1. 

Next we consider the regime |a| > 1, then X{a) < ci|a|“^/2 and we have 


m 


I < 


+ 


l(|5-g| > l/2)ds 
(s)3/2(|g| — y^)|a — s|2-“ 
f l(|a — s| < l/2)ds 


(s)V2(|g| _ [(5 _ s)2 + (A2|5|-V2 + r,)2] 


The first integral is bounded by Ca{a) — ^J2\a\) by using that (a) ~ (a). The 

second integral is bounded by 

_1_ _d£_^ ^^a-2(i-'^) 

{ay/mq\-V^) J-l/2 [s 2 + (^2|5|-1/2 + ^)2]W2 “ (a)a/2(|g|_^) 

and this completes the proof of (IB.21) . □ 
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1 


We now prove the more accurate estimate (im . We have 
A2 A2 


Im 


|a — cu(p) — \‘^Z{e{p)) + t] a — e{p) — \‘^Tl{e{p)) — i{\‘^Z{e{p)) + rf) 

From the resolvent identity and with the notations e = e(p), a = a — \^TZ{a), the last 
term equals to (/) + (//) + {HI) with 

A2 . 1 


^ ^ XH{a) + p a — e — i{XH{a) + p) 
A2 A2(J(e)-J(5))^ 


(B.7) 


{II) := - 


{III) := - 


XH{a) + p XH{e) + p 


-Im 


1 


A2 


XH{a) + p 


Im 


1 


a — e — A^0(e) — ip 

A2(0(a) -0(e)) 


a — e — i{XH{a) + p) a — e — A20(e) — ip 


Our goal is to estimate f \B{p — ( 7 )p[(/) + {II) + (///)] dp. 

We recall two continuity properties of 0£(a,r) from Lemma 3.1 of m-- 

|0£(Q!,r) - 0e(Q!,r')| < c\ |r| - \r'\ \ (B.8) 

|0e(a, r) — 0e/(a^,r)| <C{\£ — £'\ + \a — (B.9) 


if e > e' > 0. In estimating the integral of term (I), we first use (IB. 81) to change q to q 
with e{q) = a, then we use (IB. 91) with e' —> 0 + 0 and e = A^X(5) + p = 0{X^) to obtain 


\B{p-q)\‘^{I) dp = 


Z{a) + 0{X) + 0{\a-e{q)\^/^) 


XH{a) + p 
< 1 + 0{X^p~^[X + \a- uj{q)\^/‘^]) . 


In the second term of (Ib 7I) we use (|B.5|) and we drop the positive Z{a) and X(e) terms 
in the denominators 

irrni < c(—\^ _ |A^X(e(p))+ r?| |5-e(p)|^/^ _ 

- V?]/ |5-e(p) + A2(7^(a)-7^(e(p)))|2 + |A2X(e(p))+^7|2 

< rX^(—\^ _ \a-e{p)\^/^ _ 

“ \pJ |5 — e(p) + A2(7^(a) — 77(e(p)))|2 + 7^2 ’ 

where we used that \X^{IZ{a) — 77(a))] ^ = 0{X^) <C p^ based upon (IB. 51) . 

To perform the dp integration, we distinguish two regimes depending on whether 
|a —e(p)| is bigger or smaller than KX^ for a sufficiently large fixed K. When |a —e(p)| > 
KX^, then A2|77(a) — 77(e(p))| < ||a — e(p)|, hence 


i(,/)i<cv(X)" 


,-1/2 


p J \a — e(p)| + p ’ 
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and, by using a bound analogous to (ira . the corresponding integral is bounded by 


caV'/' [ < OiX%-^/^\\ogv\) . 

J \a-e{p)\+v 

When \a — e{p)\ < KX'^, then we can trivially estimate |(//)| < C{X‘^p~^)^ and after 
the co-area formula, the volume factor is given by 

1(|5 - s| < KA^)V^5(s)ds = O(A^) , 

with 

S{e) := [ \B{V^{(l)r - (l)))\‘^d(l) 

where (jir € 5''^“^ is fixed. Recalling the properties of S'(e) from the proof of Lemma 
3.2 in cni, we see that the contribution to the integral f \B{p — q)\‘^{II) dp is of order 

0((AV^)^)- 

Finally, the last term in (jEZl) is estimated as 



|(m)| <CA' 


a 


-e|V2 


r] J \d — e\ + rj \a — e + X^'R.{e)\ -|- p 


In the regime where \a — e{p)\ > KX? (with some large K) we obtain 




1 


p J \a — e\+ p 


and after integration we collect 0{X'^p~‘^/‘^\ logT/|). In the regime where |a —e(p)| < K)? 
we have < 0{X?p~^) and the volume factor is O(A^), therefore the integral is 

C>(A^??“^). Collecting the error terms we arrive at the proof of Lemma I2.1L □ 


B.2 Proof of Lemma 16.11 

We can assume that / is a real function and write f{p) = {p)~‘^'^g{p) with ||p|| 2 do < 
C\\f\\M,l. We can restrict the integration regime in (jns to IpI < A ^ since the con¬ 
tribution of the outside regime is 0(A^'^) by a Schwarz inequality (to separate the two 
denominators) and a trivial application of Lemma fB.ll with a = 0. This large momentum 
cutoff will be done with the insertion of a function xiMp)) with a smooth, compactly 
supported Xi X = 1 on [—1,1]. 

We can also assume that la —e(p—r)| < A, \(5 — e{p+r)\ < A, otherwise at least one of 
the denominator can be estimated by 0(A“^) and the other one integrated out by (I2.hjl to 
give 0(A| log A|). Since |a —e(p —r)| > |a —e(p)| —Cdpj-l-IrDIrl > |a —e(p)| — 
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we obtain that \a — e{p)\ < 2A and similarly \f5 — e{p)\ < 2A, in particular \a — j3\ < 4A 
and I 7 — e{p)\ < 2A. 

We replace the hrst denominator of (EH) hy a — e{p)+p-r — — ip. The error 

term of this replacement, by the resolvent expansion, is bounded by 



_ A^/^x(A(p))|/(p)|dp 

— r) — ip\ \a — e{p) + p ■ r — A^ 0 ( 7 ) 


ip\ \f3 — uj{p + r) + ip\ ’ 


(B.IO) 


where we have used the estimate 


a 


— e{p — r) — \^6{p — r) — ip = a — e{p) + p- r — X^Q{y) — ip + 0{X^^‘^) , 


which follows from the above restrictions on the integration domain |r| < and 

the Holder continuity (IB.51) . To estimate the error term (|B.10|1 . we bound the (3 denom¬ 
inator trivially by p~^ and use the Schwarz inequality to separate the remaining two 
denominators 


1 


< 


1 


+ 


1 


a — ijj{p — r) — ip\ \a — e{p) | \a — u;{p — r) — ip\‘^ \a — e{p) . P 


The integral of the first term can be bounded by Lemma |B3 with a = 0; for the second 
term we rewrite e{p) + p ■ r = e{p + r) — and use (12. 8 j) after a shift in a and p. We 
arrive at H with 


Hn — 


/ 


A^x(A(p))/(p)dp 


a — e{p) + p ■ r — X‘^Q{y) — ip)\P — e{p) — p ■ r — X‘^Q{'y) + ip 


To compute XIq we can choose a coordinate system where the vector r points in the n-th 
direction: r = |r|(0,..., 0,1). We can write 


Hn — A' 


I 

Jr 


dp|| 


/ 

Jr-' 


a - P + 2\r\ p|| - 2i[X^I{'y) + p] 

1 


-P — e{p) — p|| |r| — A^ 0 ( 7 ) +ip a — e{p) + py |r| — A 20 ( 7 ) — ip. 
Lemma B.2 Let F be a C^-function on R with \F{Q)\ < C{Q)~‘^ and let 

F{Q) 


(B.ll) 

x(A(j>))/(p)dp_L 


r 

Y{z) := / 
Jo 

for any z = a + ie with 0 < e < 1/2. Then 


z - Q 


dQ 


\Y{z) -Y{zp\ < C'|F(0)| llogz-logz'l -I-|z - z'||loge|||F|| 2 d,i , 


(B.12) 


where z' = a' + ie' and e > e' > 0. 















Proof. This lemma is essentially Lemma 3.10 in For completeness we recall the proof. 
Choose a branch of the complex logarithm on the upper half plane and use integration 
by parts: 


\Y{z)-Y{z')\ < |F(0)||logz 


poo 

\ogz'\+ / F'{Q)[\og{z-Q)-\og{z' 

JO 


Q)]dQ 


The second term is estimated by 



where r(z, z') is any path in the upper half plane that connects z and z' and d|^| is the 
arclength measure. A simple exercise shows 


log(Im^)| . 

Choose a path from z = a Y ie to a' + is then to a' + is' along straight line segments. 
After integration we obtain (IB.1211 . □ 


\F'm 

ie-gi 


dQ < C\\F\\2d,i 


We now change the denominators in the square bracket in (IBTTI) to 7 — e{p) ± 
ir]. This requires a change of order 0(A) in the denominators using the estimates on 
\a — 7 I and P|||r|. With the help of Lemma E21 such change yields an error of order 
A^l log A| 7 “^||/|| 2 d,i in flo. After these changes, we can remove the cutoff x(A(p)) at a 
price of as before, and we have 


f M\ f 

\ 1 1 ] 

/r a - /3 + 2|r|p|| - 2i[X‘^I{-y) + 7 ] J^d-i 

-7 — e{p) + irj 7 — e{p) — ip. 


flo = A^ 


modulo negligible errors involving ||j|| 4 d,i- 
The inner integral is evaluated as 


f{p)dpi. 

(B.13) 


I := i Im 


f*(u,p\\)u 2 du /" 1/o 

/ - 1 2 1 — f{u,P\\)-= f{u^^^ 6 ,p\i)d 9 . 

Jo 7 - 2 l^j| - + W 


For p|| in the range 7 — < —A, we have 


1^1 < 


g^||/||4d,0 

Jo |A + inp 


0(7A-1) , 


by using the decay of f* inherited from /. The contribution of this regime to Oq is 
therefore of order A| log A|, and hence negligible. In the regime I 7 — \p\\\ < A one can 
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estimate |/| < C{p\\) The dp||-volume is at most so the contribution of this 

regime to is at most of order and hence also negligible. 

Finally, we can concentrate on the regime 7 — > A. We can use the estimate (for 

e > e' > 0) 

Im f ^ dx = — 7 rff( 0 ) + 0{e je) + 0(e'\ log£^|) 

J _ £ X “h XS 

if 5 G with a bounded derivative. We obtain, with 6 = 7 — , e' = rj, that 

I = -27ri(27 -p|) V/*(27 -p|,p||) + 0(A^+^'‘) , 

where the error is integrable in py. Therefore it is negligible in Qq. Substituting the 
main term into (IB. 131) . we obtain the main term in (lOl) . □. 

C General estimates on circle graphs 

We define four operations on a partition given on the vertex set of a circle graph on 
N vertices and we estimate how the Fl-value of the partition changes. Operation I was 
already defined in Section 9 of uni, here we repeat the definition and the corresponding 
estimate for convenience. 

Operation I: Breaking np Inmps 

Consider a Feynman graph on N vertices (Section EU). Given a partition of the set 
V\{0,0*}, P = {P^ : p G .^(P)} £ PV) we define a new partition P* by breaking up one 
of the lumps into two smaller nonempty lumps. Let = Py' U Pyu with Pyi D Pyu = 0 
and P* = {Py', Py", Pfj, : € I(P) \ particular /(P*) = /(P) U {u' , u"} \ {i'} 

and m(P*) = m(P) + 1. The following estimate was proven in Lemma 9.5 of m- 

Lemma C.l With the notation above, we have 

E{^)g{P,u,a) < / dr Ef^^jg(P*,u*{r,u),OL) , 

J\r\<N( 

where the new set of momenta u* = u*(r, z^) is given by u* := u^, p G L(P) \ {n} and 
Uyi = Uy — r, u*„ = r. In our estimates we will always have N < 2K and then 

supFlp)g(P,u, q) < AFlp)gSup(P*,u,«) 

U U 

with A := [CK(Y = (see (12.181) and (14. 8j) ). 
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Operation II: Removing the lump of a single vertex 

Let V G V \ {0,0*} be a vertex and let P e Pv such that Pa = {uj for some 
a G I{P), i.e. the single element set {uj is a lump. Define V* := V \ (uj, £(V*) := 
C{V) U {(u — 1, u + 1)} \ {(u — 1, v), {v, V + 1)1, i.e. we simply remove the vertex v from 
the circle graph and connect the vertices v — l,v + 1. Let P* G Vv*, P* := P \ { {u} } 
be P after simply removing the lump {uj. In particular, /(P*) = /(P) \ {cr}. 


Lemma C.2 With the notations above 


supL^(*)g(P,u,Q:) < CAr? ^ sup.E(^)g+i(P*, u*, a) . (C.l) 

u u* 

If both neighbors of 0*, v v', form single lumps in P, then both of these lumps can 
be simultaneously removed to obtain a partition P* := P \ {{t'j, with the estimate 

supi?*g(P,u, a) < CX^ supLig+ 2 (P*, u*, q:) . (C.2) 

u u* 


Proof. The factor A in the estimate is due to the fact that each vertex (apart 

from 0 and 0*) carries a factor A and |V*| = |V| — 1. Let Pi^ be the lump of the vertex 
V — 1 right before to v in the circular ordering and assume u — 1 7 ^ 0, 0* (otherwise we 
consider v + 1 and the proof is slightly modified). We use the trivial bound 


|ae„- - U}{we,_) + ir]\ 


< rj 


-1 


(C.3) 


(recall that Im to < 0 from Lemma 3.2 of ^U]) and the bound 


\BiWe,+ - We^_)\ “ W, 


< 


C 




(C.4) 


(uniformly in tCe„_) to obtain the necessary decay between the two newly consecutive 
momenta. The same bound holds if some of the B{-) on the left hand side is replaced 
with (•)“^'^ due to the set Q. Now we integrate to obtain a new delta function from 


dhiwe„-)s(we 


- We 


- Ur. 


+ Z] 

e&L±(Pu) : e7^e„ 


±We - U, 


and clearly 


< ^ (we.,+ + Z - (Ua + Uufj 

eeL±{Pir) ■■ e^e„_ 

We,+ + Z = Z • 

e£L±{Pe) ■■ e^ev- e£L±{P*) 
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The new auxiliary momentum associated to Py is and disappeared, so the 

sum of the auxiliary momenta remain unchanged, and (ESI) continues to hold. This 
proves dnu). 

For the proof of (irr2i) . if V and v' are the vertices on both sides of 0*, then they 
can be removed and their neighbours can be connected directly to 0* yielding a non- 
truncated value of a graph with two vertices less. This gives a factor A^. The appropriate 
redefinition of the auxiliary momenta is straightforward. □ 

Operation III: Removing half of a gate 

Let v,v + 1 G V \ {0,0*} two subsequent vertices and let P G "Py such that v = 
V + 1 (mod P). In the main application this will arise when v,v P 1 are connected and 
form a gate. Define V* := V\{u+1}, TCF*) := £(V)U{(u, u+2)}\{(u, u+1), (u+1, u+2)}, 
i.e. we simply remove the vertex v + 1 from the circle graph with the adjacent edges and 
add a new edge between the vertices v,v + 2. Let P* G Py* be identical to the partition 
P except that u + 1 is simply removed from its lump. In particular, /(P) = /(P*). 

Lemma C.3 With the notations above 

E{^)g(P,M,cx) <CX\logr]\ P(*)g+i(P*,n,Q) . 

Proof. Note that the momentum of the edge between v and u+1 does not appear 
in the delta functions in the definition of P(^,)g(P, u, a) (see (01) 1. Before integrating 
out this momentum in ESI), we use the bound 

\B{We,+ -We^_)\ ) | 

{We^+ - {We,_ - 

to ensure the decay between the momenta We^_ and tCe(„+ 2 )_j that are consecutive in the 

new graph. The same bound holds if some of the B{-) is already replaced with 
The integration of rce„_ yields C\ logr/| by using (12.bl) . □ 

Operation IV: Removing an isolated gate 

Let v,v + 1 G V \ {0,0*1 be two subsequent vertices and let a partition P G Py such 
that V = v + 1 (mod P). Define V* := V \ {u, v + 1}, £(V*) ;= P(V) U {(u — 1, u + 2)} \ 
{(u — 1, u), (u, u + 1), (u + 1, u + 2)1, i.e. we simply remove the gate. Let P* G Py* be P 
after removing the lump {v,v + 1}. Combining Operations III and II, we immediately 
obtain: 
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Lemma C.4 With the notations above 


sup-E(^)g(P,u,Q!) < CA^r/ log77| sup£’(*)g+2(P*,u*,a) . □ 

U U* 

Note that this bound is not optimal. The removal of a gate affects the value of the graph 
only a by constant factor, but the corresponding estimate is more complicated and we 
do not aim at optimizing the value of k. 
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